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COGNITION ANALYSIS 



CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] This application claims priority to U.S. Application Serial 

No. 60/492,053, ifiled on 1 August 2003, the contents of which is hereby incorporated 

by reference in its entirety. 

10 RESERVATION OF COPYRIGHTS 

[0002] The disclosure of this patent document contains material which is 

subject to copyright protection. The copyright owner has no objection to the facsimile 
reproduction by anyone of the patent document or the patent disclosure, as it appears 
in the Patent and Trademark Office patent file or records, but otherwise reserves all 
1 5 copyright rights whatsoever. 

BACKGROUND 

[0003] Neuroimaging has been used to detecting abnormalities in individuals 

that suffer from neuropsychiatric disorders. However, the conventional methods for 
evaluating neuropsychiatric disorders rely on outwards signs or "exophenotypes" of 
20 iUness. The American Psychiatric Association's Diagnostic Statistical Manual (DSM- 
IV, 1994) is an example of a diagnostic procedure that uses such exophenotypes. 
Further, a number of psychiatric disorders may be caused, at least in part, by genetic 
components, the vast majority of which remain unidentified. 

SUMMARY 

25 [0004] Methods for evaluating information about the structure and fimction of 

neural circuits in the brain can be used for diagnosis and gene identification and, 
accordingly, are of particular importance for medicine, pharmacology, and society. 
Many of the methods and data management features described herem consoUdate 
relationships within multi-dimensional complex data sets, e.g., data sets that include 

30 systems biology measures, such as those obtained from neuroimaging, and, optionally 
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5 also genetic measures, e.g., from the same individuals. Generally, this process can be 
applied to any multi-variable space, not just to quantitative measures from the brain or 
genome, hi the context of neuroimaging and genetics, this process is one that has 
direct implications for the identification of genes for susceptibihty and/or resistance to 
fimctional brain illness, e.g., a behavioral or cognitive illness. 

10 [0005] Accordingly, in one aspect, the invention features a datastnxcture that 

includes: a) genetic information that describes a plurality of genetic markers of a 
subject or a reference to such information; and b) a systems biology map of the 
subject or a reference to such a map, e.g., wherein the map includes information about 
neural circuit function in the brain. The datastructure can be encoded in machine- 

15 accessible media or memory. The datastructure can also be transmitted, e.g., as a 
signal (e.g., a modulated or encoded) or other conununication, e.g., electronically or 
digitally. 

[0006] La one embodiment, the systems biology map includes structural 

information (e.g., only structural information). 

20 [0007] In one embodiment, the systems biology map includes fimctional 

information. For example, the systems biology map includes information about 
activity in a plurality of brain regions in at least one mental process, e.g., a paradigm, 
e.g., in at least two, three, four, or five paradigms. The paradigm, typically, includes 
an extemal framework with which a mental process interacts. For example, the 

25 mental process can be made to interact with an extemal stimulus, an extemal request, 
an extemal task, or an extemal sequence. 

[0008] In one embodiment, the pluraUty of brain regions includes at least five, 

ten, twenty, thirty, forty, fifty, or sixty brain regions. For example, at least one, ten, 
twenty, or thirty of the brain regions of the pluraUty are selected from Table 1. 
30 Regions can be defined by structural and/or fimctional features. Subregions or 
smaller volumes than the exemplary regions in Table 1 can also be used, as can 
regions that are defined by larger volumes and that encompass one or more of the 
exemplary regions. 
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5 [0009] In one embodiment, the information for each of the brain regions is 

independent of reference to a coordinate frame. For example, the brain regions can be 
identified by a numerical index (e.g., an index values for each of a set of predefined 
regions) or by text (e.g., a categorical reference) or an indirect reference (e.g., use of 
pointers and hyperlinks). In another embodiment, one or more the brain regions can 
10 be identified by reference to a coordinate frame, e.g., Talairach coordinates. For 

example, the information is not indexed voxel by voxel so as not to be in a form of a 
raster, e.g., the information is non-rasterized. 

[001 0] In one embodiment, the paradigm interacts with the informational 

backbone for motivation, e.g., it evokes at least one region in the informational 

IS backbone for motivation. In one embodiment, the paradigm interacts with 

mechanisms for representation and convergence, feature evaluation, probability 
assessment, outcome processing, valuation, reward/aversion processing, 
coimterfactual comparisons, and memory. In one embodiment, the paradigm interacts 
with mechanisms for selection of objectives for fitness, mechanisms for selection of 

20 behavior, or information processing (e.g., reception). In one embodiment, the 
paradigm interacts with mechanisms for language and symbol processing, 
mechanisms for commimication, and/or mechanisms for social behavior. 

[001 1] In one embodiment, the systems biology map can include information 

obtained by imaging, e.g., neuroimaging, e.g., tomography, e.g., MRI, fMRI, MEG, 
25 fCT, OI, SPECT, or PET system. Neuroimaging can including imaging at least one 
region of the brain or central nervous system. 

[0012] In an embodiment in which there is information for least two 

paradigms, these paradigms may interact with overlapping, but non-coextensive 
regions of the brain, e.g., each paradigm may interact with at least one region that is 
30 not activated in another paradigm by a normal subject. Exemplary paradigms include: 
a social reward paradigm, a CPT / probability paradigm, a physiological aversion / 
pain paradigm, a mental rotation paradigm, an emotional faces paradigm, and a 
monetary reward paradigm. Other paradigms can also be used. For example, another 
paradigm which interrogates the informational backbone for motivation or other areas 
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5 described herein, e.g., an area interrogated by one of the above paradigms can be 
used. 

[0013] In one embodiment, the information about activity for at least one of 

the regions includes deviations from a reference (e.g., percentage difTerences, ratios, 
and subtractive values). 

10 [0014] In one embodiment, the systems biology map includes one or a 

plurality of matrices, each matrix including information about neural activity in a 
plurality of defined brain regions during different paradigms. In another embodiment, 
the map includes a similar or identical set of information, but is stored or represented 
in another form, e.g., as text, graphic, e.g., as a vector, table, etc. In one embodiment, 

15 two matrices are used, one including information about individual activation regions 
and another about population activation regions. 

[0015] For example, the plurality of genetic markers includes markers on at 

least two, three, four, five, six, ten, twelve, or fifteen different, non-homologous 
chromosomes. In one embodiment, theplurality of genetic markers includes markers 
20 on each autosome, e.g., at least one, two, five, ten, twenty, or fifty markers on each 
autosome. For example, at least 20, 50, or 70% of the markers can be spaced closer 
than 500, 50, 20, 10, or 2 Mb to another marker c» 200, 100, 50, 20, or 10 cM to 
another marker. 

[0016] The genetic information can includes information about, e.g., 

25 nucleotide identity for a plurality of genetic markers, methylation status for a plurahty 
of genetic markers, parental origin for one or a plurality of genetic markers, chromatin 
structure or accessibility for one or a plurality of genetic markers, a haplotype, 
microsatellite marker, sequence tagged site, SNP, a chromosomal deletion, inversion, 
transversion, rearrangement, trisomy, or other chromosomal abnormality. 

30 [0017] In another aspect, the invention features a datastructure that includes: a 

systems biology map of a subject wherein the map includes quantitative information 
about neural circuit function in the brain. For example, the information indicates 
function of a plurality of regions of the brain during a plurality of mental processes. 
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5 [001 8] In one embodiment, the systems biology map includes structural 

information (e.g., only structural information) 

[0019] In one embodiment, the systems biology map includes functional 

information. For example, the systems biology map includes information about 
activity in a plurality of brain regions in at least one mental process, e.g., a paradigm, 

10 e.g., in at least two, three, four, or five paradigms. In one embodiment, the plurality 
of br?dn regions includes at least five, ten, twenty, thirty, forty, fifty, or sixty brain 
regions. For example, at least one, ten, twenty, or thirty of the brain regions of the 
plurality are selected firom Table 1 . Subregions or smaller volumes than the 
exemplary regions in Table 1 can also be used, as can regions that are defined by 

1 5 larger volumes and encompass one or more of the exemplary regions. 

[0020] In one embodiment, the information for each of the brain regions is 

independent of reference to a coordinate firame. For example, the brain regions can be 
identified by a numerical index (e.g., an index values for each of a set of predefined 
' regions) or by text (e.g., a categorical reference) or an indirect reference (e.g., use of 
20 pointers and hyperlinks). In another embodiment, one or more the brain regions can 
be identified by reference to a coordinate fi^ame, e.g., Talairach coordinates. For 
example, however, the information is not indexed voxel by voxel so as not to be in a 
form of a raster, i.e., the information is non-rasterized. 

[0021] In one embodiment, the paradigm interacts with the informational 

25 backbone for motivation, e.g., it evokes at least one region in the informational 
backbone for motivation. In one embodiment, the paradigm evokes at least one 
region in the informational backbone for motivation. In one embodiment, the 
paradigm interacts with mechanisms for representation and convergence, feature 
evaluation, probability assessment, outcome processing, valuation, reward/aversion 
30 processing, counterfactual comparisons, and memory. In one embodiment, the 
paradigm interacts with mechanisufis for selection of objectives for fitness, 
mechanisms for selection of behavior, or information processing (e.g., reception). In 
one embodiment, the paradigm interacts with mechanisms for language and symbol 
processing, mechanisms for communication, and/or mechanisms for social behavior. 
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5 [0022] The systems biology map can include infonnation obtained by 

imaging, e.g., neuroimaging, e.g., tomography, e.g., MRI, fMRI, MEG, fCT, OI, 
SPECT, or PET system. 

[0023] In an embodiment in which there is infonnation for least two 

paradigms, these paradigms may interact with overlapping, but non-coextensive 
10 regions of the brain, e.g., each paradigm may interact with at least one region that is 
not activated in another paradigm by a normal subject. 

[0024] Exemplary paradigms include: a social reward paradigm, a CPT / 

probabihty paradigm, a physiological aversion / pain paradigm, a mental rotation 
paradigm, an emotional faces paradigm, and a monetary reward paradigm. Other 
15 paradigms can also be used. For example, another paradigm which interrogates the 
informational backbone for motivation or other areas described herein, e.g., an area 
interrogated by one of the above paradigms can be used. 

[0025] In one embodiment, the mfbrmation about activity for at least one of 

the regions includes deviations Scorn a reference (e.g., percentage differences, ratios, 
20 and subtractive values). 

[0026] In one embodiment, the systems biology map includes a plurality of 

matrices, each matrix including information about neural activity in a plurality of 
defined brain regions during different paradigms. In another embodiment, the map 
includes a similar or identical set of information, but is stored or represented in 
25 another form, e.g., as text, graphic, e.g., as a vector, table, etc. 

[0027] The datastructure can further include genetic information that describes 

a plurality of genetic markers of the subject or a reference to such information. For 
example, the plurality of genetic markers includes markers on at least two, three, four, 
five, six, ten, twelve, or fifteen different, non-homologous chromosomes. In one 
30 embodiment, the plurality of genetic markers includes markers on each autosome, 
e.g., at least one, two, five, ten, twenty, or fifty markers on each autosome. For 
example, at least 20, 50, or 70% of the markers can be spaced closer than 500, 50, 20, 
10, or 2 Mb to another marker or 200, 100, 50, 20, or 10 cM to another marker. 
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5 [0028] The genetic information can includes information about, e.g., 

nucleotide identity for a plurality of genetic markers, methylation status for a plurality 
of genetic markers, parental origin for one or a plurality of genetic markers, chronaatin 
structure or accessibility for one or a pliurality of genetic markers, a haplotype, 
microsatellite marker, sequence tagged site, SNP, a chromosomal deletion, inversion, 
10 transversion, rearrangement, trisomy, or other chromosomal abnormality. 

[0029] The datastructure can be encoded in machine-accessible media or 

memory. The datastructure can also be transmitted, e.g., as a signal (e.g., a 
modulated or encoded) or other conmiunication, e.g., electronically or digitally. 

[0030] In another aspect, the invention features a datastructure that includes: a 

15 systems biology map of a subject wherein the map includes a plurality of vedues 
corresponding to a set of continuous variables, wherein the variables of the set 
correspond to different regions of the brain, and the values that correspond to the 
variables indicate fimction of respective regions during a mental process. 

[0031] In one embodiment, the systems biology map includes structural 

20 information (e.g., only structural information) 

[0032] In one embodiment, the systems biology map includes functional 

information. For example, the systems biology map includes information about 
activity in a pluraUty of brain regions in at least one mental process, e.g., a paradigm, 
e.g., in at least two, three, four, or five paradigms. In one embodiment, the plurality 
25 of brain regions includes at least five, ten, twenty, thirty, forty, fifty, or sixty brain 
regions. For example, at least one, ten, twenty, or thirty of the brain regions of the 
plurality are selected firom Table 1. Subregions or smaller volumes than the 
exemplary regions in Table 1 can also be used, as can regions that are defined by 
larger volumes and encompass one or more of the exemplary regions. 

30 [0033] In one embodiment, the information for each of the brain regions is 

independent of reference to a coordinate firame. For example, the brain regions can be 
identified by a numerical index (e.g., an index values for each of a set of predefined 
regions) or by text (e.g., a categorical reference) or an indirect reference (e.g., use of 
pointers and hyperlinks). In another embodiment, one or more the brain regions can 
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5 be identified by reference to a coordinate frame, e.g., Talairach coordinates. For 
example, however, the information is not indexed voxel by voxel so as not to be in a 
form of a raster, i.e., the information is non-rasterized. 

[0034] In one embodiment, the paradigm interacts with the informational 

backbone for motivation, e.g., it evokes at least one region in the informational 

10 backbone for motivation. In one embodiment, the paradigm evokes at least one 
region in the informational backbone for motivation. In one embodiment, the 
paradigm interacts with mechanisms for representation and convergence, feature 
evaluation, probability assessment, outcome processing, valuation, reward/aversion 
processing, counterfactual comparisons, and memory. In one embodiment, the 

15 paradigm interacts with mechanisms for selection of objectives for fitness, 

mechanisms for selection of behavior, or information processing (e.g., reception). In 
one embodiment, the paradigm interacts with mechanisms for language and symbol 
processing, mechanisms for communication, and/or mechanisms for social behavior. 

[0035] In one embodiment, the systems biology map is condensed relative to a 

20 native dataset (e.g., arasterized dataset), e.g., at least 10, 10^ 10^, 10\ 10^ or 10^ 
fold. 

[0036] In an embodiment in which there is information for least two 

paradigms, these paradigms may interact with overlapping, but non-coextensive 
regions of the brain, e.g., each paradigm may interact with at least one region that is 
25 not activated in another paradigm by a normal subject. 

[0037] Exemplary paradigms include: a social reward paradigm, a CPT / 

probability paradigm, a physiological aversion / pain paradigm, a mental rotation 
paradigm, an emotional faces paradigm, and a monetary reward paradigm. Other 
paradigms can also be used. For example, another paradigm which interrogates the 
30 informational backbone for motivation or other areas described herein, e.g., an area 
interrogated by one of the above paradigms can be used. 

[0038] In one embodiment, the information about activity for at least one of 

the regions includes deviations firom a reference (e.g., percentage differences, ratios, 
and subtractive values). 



-8" 



wo 2005/020788 



PCT/US2004/025146 



5 [0039] In one embodiment, the systems biology map includes a pliirality of 

matrices, each matrix including information about neural activity in a plurality of 
defined brain regions during different paradigms. In another embodiment, the map 
includes a similar or identical set of information, but is stored or represented in 
another form, e.g., as text, graphic, e.g., as a vector, table, etc. 

10 [0040] The datastructure can further include genetic information that describes 

a plurality of genetic markers of the subject or a reference to such information. For 
example, the plurality of genetic markers includes markers on at least two, three, four, 
five, six, ten, twelve, or fifteen different, non-homologous chromosomes. In one 
embodiment, the plurality of genetic markers includes markers on each autosome, 

15 e.g., at least one, two, five, ten, twenty, or fifty markers on each autosome. For 

example, at least 20, SO, or 70% of the markers can be spaced closer than 500, 50, 20, 
10, or 2 Mb to another marker or 200, 100, 50, 20, or 10 cM to another marker. 

[0041] The genetic information can includes infoiraation about, e.g., 

nucleotide identity for a plurality of genetic markers, methylation status for a plurality 
20 of genetic markers, parental origin for one or a plurality of genetic markers, chromatin 
structure or accessibility for one or a plurality of genetic markers, a haplotype, 
microsatellite marker, sequence tagged site, SNP, a chromosomal deletion, inversion, 
transversion, rearrangement, trisomy, or other chromosomal abnormality. 

[0042] In one embodiment, the datastructure further includes c) information 

25 that is an index corresponding to the subject. For example, the index can be 
randomized, encrypted, or anonymous. In another example, the index directly 
identifies the subject (e.g., name, social security number etc). In one embodiment, the 
index associates the subject with famihal or other pedigree information. 

[0043] The datastructure can be encoded in machine-accessible media or 

30 memory. The datastructure can also be transmitted, e.g., as a signal (e.g., a 

modulated or encoded) or other communication, e.g., electronically or digitally. 

[0044] The invention also features database including: a plurality of records, 

wherein each record of the plurality includes a datastructure described herein or other 
datastructure which is condensed relative to native information (e.g., rasterized data) 
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5 obtained from subjects at a plurality of time points. In one embodiment, the 

datastructure is accessible to statistical analysis (e.g., uncompressed) and enables 
phenotypic classification of subjects. 

[0045] In one embodiment, the records of the plurality include records for a 

plurahty of unrelated individuals and records for at least one biological family 
10 member of each of the plurality of unrelated individuals. For example, at least 5, 10, 
20, 30, or 50% of the database can include records for individuals for which there is 
.also a record for a biologically related family member. 

[0046] In one embodiment, the database includes records for at least 50, 100, 

200, 500, 1000, 3000 or 30,000 human subjects, or ranges therebetween. In one 
15 embodiment, the database includes records for individuals from different populations 
(e.g., ethnic populations, e.g., at least two, three, or four different continents, e.g., 
Caucasians, Africans, Polynesians, Native Americans, and Asians). 

[0047] In one embodiment, the database includes records for at least 50, 100, 

200, 500, 1000, 3000 or 10,000 human subjects who each have a clinical diagnosis of 

20 a neurological and/or psychiatric disorder, e.g., schizophrenia, manic dqpression, 
bipolar disorder, addictions (e.g., substance abuse, gambling, etc.), obsessive- 
compulsive disorder, anxiety/paranoia, autism, schizo-affective disorder, delusional 
disorder, psychosis, antisocial personaUty disorder, or anorexia/bulimia nervosa. For 
example, the database can include at least 50, 100, 200, 500, 1000, 3000 or 10,000 for 

25 a single disorder. 

[0048] For example, the datastructure is a condensed form of a native dataset, 

e.g., (arasterized dataset). For example, the datastructure is condensed at least 10, 
10^, 10\ 10^ 10^ or 10^ fold. 

[0049] In one embodiment, the datastructure further includes genetic 

30 information or a reference to such information. The datastructure can include other 
features described herein. 

[0050] In another aspect, the invention features a method that includes: 

providing a database that includes information about brain activity (e.g., structural 
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5 and/or functional information) for each of a plurality of subjects (e.g., a database 
described herein); and 

[005 1 ] classifying the subjects based on the information. 

[0052] In one embodiment, the classifying includes selecting a subset of 

variables, and sorting the subjects as a function of the variables of the subset. For 
10 example, the subset of variables can be selected based on the information content 
(e.g., relative information content) of each of the variables. For example, the subset 
of variables can be selected based on correlations (e.g., autocorrelations) among the 
variables. In one embodiment, each variable is associated with an activity of a brain 
region and a mental process, e.g., a paradigm. 

15 [0053] In one embodiment, the classifying includes generating, evaluating, or 

characterizing a tree, e.g., a binary tree. For example, each node of the tree 
corresponds to a variable associated with a particular region of the brain and a mental 
process. 

[0054] In one embodiment, the classifying is recursive. 

20 [005 5] In one embodiment, the plurality of subj ects includes at least 50, 1 00, 

200, 500, 1000, or 3000 human subjects. 

[0056] In one embodiment, the classifying includes an association rule . 

algorithm. For example, the association rule algorithm is non-parametric. 

[0057] In one embodiment, the classifying includes a classification tree 

25 analysis, hierarchical clustering, Bayesian clustering, k-means clustering, self- 
organizing maps, and/or shortest path analysis. 

[0058] In one embodiment, the method further includes comparing genetic 

information among subjects of at least one class, e.g., evaluating a statistic for 
association of one or more genetic markers among the subjects of the at least one 
30 class. 

[0059] In one embodiment, the information includes quantitative volumetric 

data evaluated by tomography, e.g., MRI, e.g., fMRI or mMRI. The quantitative 
volumetric data can include a plurality of matrices. 
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5 [0060] In one embodiment, the subjects are social non-human animals, e.g., 

non-human primates or voles. In one embodiment, the subjects are voles. 

[0061] In another aspect, the invention features a method that includes: 

providing a database that includes quantitative information about brain function for 
each of a plurality of subjects; and identifying, e.g., objectively identifying, a subset 
10 of subjects from the plurality of subjects according to similarity of brain function. In 
one embodiment, a plurality of subsets are objectively identified. 

[0062] In one embodiment, the identifying includes selecting, e.g., objectively 

selecting, a subset of quantitative variables whose values vary among the pluraUty of 
subjects. 

15 [0063] In one embodiment, the method furflier includes receiving additional 

quantitative information about brain function for at least one additional subject, and 
evaluating v^^hether the additional subject is a member of the identified subset. 

[0064] For example, the identifying includes generating one or more 

association rules that model the subset; a decision tree that models the subset; and a 
20 probability function that models the subset. In one embodiment, the database includes 
systems biology maps. For example, the systems biology maps includes values 
determined evaluating subjects during at least two different mental processes. 

[0065] In another aspect, the invention features a data-tree that includes a 

pluraUty of nodes, wherein each non-terminal node includes (i) a reference to a 
25 variable or variable class, wherein the variable or variable class is a parameter of 
brain function in the subject, (ii) optionally, a node level, and (iii) criterion for 
distinguish descendants of the node. 

[0066] For example, the tree is a binary tree. In one embodiment, each non- 

terminal node includes a pointer to one or more descendant nodes. 

30 [0067] In one embodiment, for at least some of the nodes of the plurality, the 

criterion is an association rule. In one embodiment, each descendant node can be 
defined by a function, e.g., a probabilistic or statistical function, that differentiates it 
from a sibling descendant node. In one embodiment, the nodes are ordered as 
function of variables that they respectively reference, e.g., as a function of 
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5 mformation content or autocorrelations for the respective variables. For example, at 
least one of the variables or variable classes refers to a brain region in a paradigm. 

[0068] In another aspect, the invention features a datastructure including a 

plurality of matrices, wherein each matrix includes functional information obtained 
during a mental process of a subject, the matrix including at least two dimensions, a 
1 0 first dimension that identifies regions of the brain, and one or more v£ilues for each 
region, wherein the values correspond to activity levels in the respective regions 
during the mental process. In one embodiment, the second dimension identifies 
lefl/right hemisphere. 

[0069] In one embodiment, the datastructure includes a first matrix that 

15 includes functional information obtained during a first paradigm and a second matrix 
that includes functional information obtained during a second paradigm. 

[0070] In one embodiment, the datastructure includes a first matrix that 

includes first values that depend on a dataset obtained by imaging the subject at 
multiple tunepoints (e.g., a native dataset, e.g., rasterized data), wherein the first 
20 values are independent of information firom other subjects, and a second matrix that 
includes second values that depend on the same dataset, wherein the second values are 
determined or are selected as a function of information firom other subjects. In one 
embodiment, the second values are selected based on location of activation centers 
detected in an aggregate of image information firom a plxurality of other subjects. 

25 [0071] In one embodiment, the first values are determined and/or selected as a 

function of location of activation centers detected by clustering signal changes firom a 
baseline, wherein the signal changes are independent of information firom any other 
subject. 

[0072] The datastructure can be encoded in machine-accessible media or 

30 memory. The datastructure can also be transmitted, e.g., as a signal (e.g., a 

modulated or encoded) or other communication, e.g., electronically or digitally. 

t 

[0073] In one aspect, the invention features a melhod that includes: providing 

(e.g., imaging or receiving) native information about brain function of a subject 
during a mental process, the information including quantitative data for signals in at 
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5 least a plurality of regions; comparing signals during the mental process to reference 
signal parameters to locate regions of activity; and populating a datastructure with 
information about signals at least in the regions of activity. The method can provide a 
systems biology map. 

[0074] Li one embodiment, the reference signal parameters is function of a 

10 baseline for the subject In another embodiment, the reference signal parameters are a 
function of signals from a population of subjects. 

[0075] In one embodiment, locating regions of activity includes clustering 

signal changes relative to the reference signal parameters. For example, the clustering 
includes defining foci in a three-dimensional coordinate space. In one embodiment, 
15 the comparing includes generating a statistical map, e.g., as a function of correlation 
between a gamma function and signal changes. The method can include other 
features described herein. 

[0076] In another aspect, the invention features a method that includes: 

providing (e.g., imaging or receiving) datasets (e.g., native or rasterized datasets) 
20 about brain function for a plurality of subjects during a mental process, the 

information including quantitative data for signals in at least a plurality of regions; 
combining information from the datasets to provide an aggregate dataset; and 
localizing regions of activity in the aggregate dataset. 

[0077] In one embodiment, the combining includes one or more of: 

25 transforming native datasets to a reference coordinate frame, averaging the native 

datasets, and producing a statistical map. In one embodiment, the localizing includes 
clustering signal changes in the aggregate dataset. 

[0078] In another aspect, the invention features a method that includes: 

providing (e.g., imaging or receiving) native datasets about brain function for a 
30 plurality of subjects during a mental process, the information including quantitative 
data for signals in at least a pliurality of regions; for each subject, producing a first 
systems biology map from the native dataset of the particular subject, wherein the first 
system biology map is independent of the native datasets from the other subjects, and 
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5 a second systems biology map that is a function of regions of activity detected in an 
aggregate dataset from the plurality of subjects. 

[0079] In another aspect, the invention features a method that includes: 

providing (e.g., imaging or receiving) information about structure and/or function of 
the brain of the subject, the information including quantitative data for at least a 

10 pliwality of regions; and objectively evaluating the information using quantitative 
criteria; and providing a diagnosis for the subject based on results of the evaluating. 
For example, the quantitative data includes information about brain function during a 
plurality of mental processes. In one embodiment, at least one mental process 
includes a paradigm, e.g., a paradigm that evokes the information backbone for 

15 motivation. An objective evaluation is typically completely free of the bias or 

potential bias of a hmnan analyst. Bias may still produced by blind or double-blind 
human analyst, because the analyst is using non-quantitative metrics to make a 
decision. 

[0080] In one embodiment, the evaluating includes comparing the information 

20 about the subject to a decision tree. For example, the comparing includes evaluating a 
probability of association for the information about the subject and one or more 
terminal nodes of the tree. In another example the comparing includes evaluating a 
probability of association for the information about the subject and each bifurcation of 
the tree. In another example, the evaluating includes evaluating a probability that the 
25 information about the subject is within a classification, wherein the classification is a 
function of quantitative activity measures for a plurality of brain regions, 

[0081] Information about the brain of the subject can be evaluated using rules 

at one or more nodes of the tree. Rules can be evaluated according to organization of 
the tree. The information about the brain of the subject can include functional 

30 information. One or more parameters (e.g., a functional or morphometric parameter) 
can be compared to boundaries of the bin. It is also possible to generate new binning 
criteria, e.g., by modifying the tree while the subject is being evaluated using the tree. 
The information evaluated by the tree can be information from a condensed 
representation of a native dataset obtained by imaging the brain. The information can 

35 be from, e.g., a matrix or a plurality of matrices. The information can be obtained 
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5 from a plurality of images and/or a plurality of paradigms. The method can include 
other features described herein. 

[0082] In another £ispect, the invention features a method that includes: 

imaging regions of the brain of a subject while at least one of the regions is active to 
obtain a native dataset (e.g., including rasterized image information) that includes 

10 information about activity in one oir more of the regions at a plurality of temporal 
instances (or receiving the native dataset); and condensing the native dataset to 
provide a condensed dataset that includes quantitative information about at least some 
of the imaged regions. In one embodiment, the condensed dataset includes 
information about one or more activity peaks in at least some of the imaged regions. 

15 In one embodiment, the condensed dataset discards time resolution for at least 10, 20, 
30, 50, 70, 80, 90, or 100% of the regions. 

[0083] In one embodiment, the regions are imaged by fMRI. In one 

embodiment, the condensing reduces data size at least 10, 10^,. 10^, 10"^, 10^ or 10^ 
fold. 

20 [0084] In one embodhnent, the condensed dataset includes information that 

can be represented as a matrix, one dimension of which differentiates among regions 
of the brain, (e.g., a region in Table 1). 

[0085] In another aspect, the invention features a method that includes: 

imaging regions of the brain of a subject during a mental process to obtain a dataset 
25 (e.g., a native dataset) that includes information about brain function; and popula1<ing 
variables in a matrix by extracting quantitative information from the dataset. In one 
embodiment, wherein the matrix includes at least two dimensions. 

[0086] In one embodiment, the first dimension resolves different regions of 

the brain. In one embodiment, the second dimension resolves the left and right 
30 hemisphere of the brain. In one embodiment, the matrix includes a third dimension. 
In one embodiment, information about one or more activations in a given region and 
hemisphere are provided at respective variables of the matrix. 

' [0087] In one embodiment, the information includes a Ust, the members of the 

Ust being stored at different positions along a third dimension of the matrix. In one 
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5 embodiment, the matrix does not provide information about time, e.g., the information 
about the one or more activations is not time-resolved. 

[0088] The imaging can include^ ^.g., neuroimaging, e.g., tomography, e.g., 

MRI, £MRI, MEG, fCT, OI, SPECT, or PET system. 

[0089] In one embodiment, the provide a systems biology map that includes 

10 functional information. For example, the systems biology map includes information 
about activity in a plurality of brain regions in at least one mental process, e.g., a 
paradigm, e.g., in at least two, three, four, or five paradigms. In one embodiment, the 
plurality of brain regions includes at least five, ten, twenty, thirty, forty, fifty, or sixty 
brain regions. For example, at least one, ten, twenty, or thirty of the brain regions of 
15 the plurality are selected fi:om Table 1 . Subregions or smaller volimies than the 
exemplary regions in Table 1 can also be used, as can regions that are defined by 
larger volumes and encompass one or more of the exemplary regions. 

[0090] In one embodiment, the information for each of the brain regions is 

independent of reference to a coordinate firame. For example, the brain regions can be 

20 identified by a numerical index (e.g., an index values for each of a set of predefined 
regions) or by text (e.g., a categorical reference) or an indirect reference (e.g., use of 
pointers and hyperlinks). In another embodiment, one or more the brain regions can 
be identified by reference to a coordinate firame, e.g., Talairach coordinates. For 
example, however, the information is not indexed voxel by voxel so as not to be in a 

25 form of a raster, i.e., the information is non-rasterized. 

[0091] In one embodiment, the paradigm interacts with the informational 

backbone for motivation, e.g., it evokes at least one region in the informational 
backbone for motivation. In one embodiment, the paradigm evokes at least one 
region in the informational backbone for motivation. In one embodiment, the 
30 paradigm interacts with mechanisms for representation and convergence, feature 
evaluation, probability assessment, outcome processing, valuation, reward/aversion 
processing, counterfactual comparisons, and memory. In one embodiment, the 
paradigm interacts with mechanisms for selection of objectives for fitness, 
mechanisms for selection of behavior, or information processing (e.g., reception). In 
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5 one embodiment, the paradigm interacts with mechanisms for language and symbol 
processing, mechanisms for communication, and/or mechanisms for social behavior, 

[0092] In an embodiment in which there is information for least two 

paradigms, these paradigms may interact with overlapping, but non-coextensive 
. regions of the brain, e.g., each paradigm may interact with at least one region that is 
1 0 not activated in another paradigm by a normal subject. 

[0093] Exemplary paradigms include: a social reward paradigm, a CPT / 

probabihty paradigm, a physiological aversion / pain paradigm, a mental rotation 
paradigm, an emotional faces paradigm, and a monetary reward paradigm. Other 
paradigms can also be used. For example, another paradigm which interrogates the 
15 informational backbone for motivation or other areas described herein, e.g., an area 
interrogated by one of the above paradigms can be used. ^ 

[0094] In one embodiment, the information about activity for at least one of 

the regions includes deviations from a reference (e.g., percentage differences, ratios, 
and subtractive values). 

20 [0095] In one embodiment, the systems biology map includes a plurality of 

matrices, each matrix including information about neural activity in a plurality of 
defined brain regions during different paradigms. In another embodiment, the map 
includes a similar or identical set of information, but is stored or represented in 
another form, e.g., as text, graphic, e.g., as a vector, table, etc. 

25 [0096] In another aspect, the invention features a method that includes: 

receiving a native dataset that includes imaged information about brain function of a 
subject; and populating variables in a matrix by extracting quantitative information 
from the native dataset. The method can be used to provide a systems biology map, 
e.g., as described herein. 

30 [0097] In another aspect, the invention features a method that includes: 

imaging regions of the brain of a plurality of subjects; transforming image 
information to a reference coordinate space; selecting a subset of regions for which 
activations are detected among the plurality of subjects; and producing a condensed 
dataset for each subject of the plurality wherein the condensed dataset is smaller than 
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5 the native dataset for each subject of the plurality and retains information about the 
selected subset of regions. In one embodiment, selecting the subset includes 
averaging the transformed image information and evaluating statistically significant 
changes relative to results of the averaging, hi one embodiment, selecting the subset 
includes selecting regions that differ from a reference (e.g., a baseline obtained prior 

10 or after the mental process). The method can include other features described herein, 
hi a related method, information is received for one or more subjects. A condensed 
dataset is produced for the subject. The condensed dataset retains information about 
one or more region in which activations are detected in the subject or in which at least 
one activation is detected among a plurahty of subjects. The retained information 

IS can be selected, e.g., as described herein. 

[0098] La another aspect, the invention features a method that includes: 

receiving functional information about neural circuit activity, the information being 
obtained by imaging a plurality of brain regions in a subject, and generating a dataset 
that associates each of a plurality of brain regions with quantitative information, 
20 wherein the quantitative information includes Usts of activation peaks (e.g., % signal 
change) and each list is associated with at least one of the brain regions, hi one 
embodiment, the hst is rank ordered, hi one example, the dataset is represented as a 
matrix. For example, members of each list are positioned or referenced in 
consecutive cells along one axis of the matrix. 

25 [0099] In another example, the dataset is represented as a vector or is stored in 

a relational database, e.g., as a table. The method can include other features described 
herein. 

[0100] M another aspect, the invention features method that includes: 

evaluating a subject to produce a first systems biology map of the subject; treating the 
30 subject; and evaluating the subject to produce a second systems biology map of the 
subject; wherein the first and second systems biology maps include quantitative 
information about brain function in a plurality of brain regions during at least one 
mental process, e.g., a paradigm. The method can be used, e.g., to evaluate a 
treatment, e.g., a candidate treatment or a previously validated treatment. 
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[0101] In one embodiment, treating the subject includes administering an 

agent to the subject. Examples of the agent include a pharmaceutical, a narcotic, an 
addictive substance, or a non-addictive substance, 

[0102] In one embodmient, treating the subject includes providing a non- 

invasive therapy to the subject. For example, the non-invasive treatment can include 
hypnosis, music, video, visual, superficial contacts, exercise, or physical pressure. 

[0103] The system biology maps can be maps described herein. For example, 

they can include information about activity in a plurality of brain regions in at least 
one mental process, e.g., a paradigm. They can include information about activity in 
a plurality of brain regions in at least two paradigms. 

[0104] In one embodiment, the systems biology map includes structural 

information (e.g., only structural information) 

[0105] In one embodiment, the systems biology map includes functional 

information. For example, the systems biology map includes information about 
activity in a plurality of brain regions in at least one mental process, e.g., a paradigm, 
e.g., in at least two, three, four, or five paradigms. In one embodiment, the plurality 
of brain regions includes at least five, ten, twenty, thirty, forty, fifty, or sixty brain 
regions. For example, at least one, ten, twenty, or thirty of the brain regions of the 
plurality are selected from Table 1 . Subregions or smaller volumes than the 
exemplary regions in Table 1 can also be used, as can regions that are defined by 
larger volimies and encompass one or more of the exemplary regions. 

[0106] In one embodiment, the information for each of the brain regions is 

independent of reference to a coordinate frame. For example, the brain regions can be 
identified by a numerical index (e.g., an index values for each of a set of predefined 
regions) or by text (e.g., a categorical reference) or an indirect reference (e.g., use of 
pointers and hyperUnks). In another embodiment, one or more the brain regions can 
be identified by reference to a coordinate frame, e.g., Talairach coordinates. For 
example, however, tlie information is not indexed voxel by voxel so as not to be in a 
form of a raster, i.e., the information is non-rasterized. 
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5 [0107] In one embodiment, the paradigm interacts with the informational 

backbone for motivation, e.g., it evokes at least one region in the informational 
backbone for motivation. In one embodiment, the paradigm evokes at least one 
region in the informational backbone for motivation. In one embodiment, the 
paradigm interacts with mechanisms for representation and convergence, feature 

10 evaluation, probability assessment, outcome processing, valuation, reward/aversion 
processing, counterfactuaJ comparisons, and memory. In one embodiment, the 
paradigm interacts with mechanisms for selection of objectives for fitness, 
mechanisms for selection of behavior, or information processing (e.g., reception). In 
one embodiment, the paradigm interacts with mechanisms for language and symbol 

15 processing, mechanisms for communication, and/or mechanisms for social behavior. 

[0108] The systems biology map can include information obtained by 

imaging, e.g., neuroimaging, e.g., tomography, e.g., MRI, £MRI, MEG, fCT, OI, 
SPECT, or PET system. 

[0109] In an embodiment in which there is information for least two 

20 paradigms, these paradigms may interact with overlapping, but non-coextensive 

regions of the brain, e.g., each paradigm may interact with at least one region that is 
not activated in another paradigm by a normal subject. 

[0110] Exemplary paradigms include: a social reward paradigm, a CPT / 

probabiUty paradigm, a physiological aversion / pain paradigm, a mental rotation 
25 paradigm, an emotional faces paradigm, and a monetary reward paradigm. Other 
paradigms can also be used. For example, another paradigm which interrogates the 
informational backbone for motivation or other areas described herein, e.g., an area 
interrogated by one of the above paradigms can be used. 

[01 11] In one embodiment, the information about activity for at least one of 

30 the regions includes deviations from a reference (e.g., percentage differences, ratios, 
and subtractive values). 

[01 12] In one embodiment, the systems biology map includes a plurality of 

matrices, each matrix including information about neural activity in a plurality of 
defined brain regions during different paradigms. In another embodiment, the map 
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S includes a similar or identical set of information, but is stored or represented in 
another form, e.g., as text, graphic, e.g., as a vector, table, etc. 

[0113] In another aspect, the invention features a method that includes: 

providing a dataset that includes quantitative information about brain activity during 
at least two paradigms for each of a plurahty of subjects; evaluating a parameter that 
10 is a continuous function of at least two components of the quantitative information, 
the at least two components being associated with different paradigms; and analyzing 
a statistic for association between the parameter and an allele for one or more genetic 
loci. For example, analyzing the statistic can include a linkage analysis, e.g., non- 
parametric linkage analysis. 

15 [01 14] 111 another aspect, the invention features methods of providing a 

population-based statistic for a brain structure. The methods include: evaluating 
images of a brain structure, e.g., for each of a plurality of subjects; aligning the 
images; and determining positional information defining an virtual brain structure 
whose structural features are based on a pre-selected probability value, the value 

20 representing the probabihty that the brain structure of one of the members of the 
plurahty is within the constraint of the virtual brain structure. For example, the 
positional information represents an isofonn surface. In one embodiment, the pre- 
selected probability value is 10, 20, 30, 40, 50, 60, 70, 80, 85, 90, 91, 92, 93, 94, 95, 
96, 97, 98, or 99%. 

25 [01 15] The plurahty of subjects can include between 5-5000 subjects, e.g., 8- 

500, 15-50, 10-35, or 51-200 subjects. For example, each subject of the plurahty has 
a common characteristic, e.g., a common behavioral trait that differs from normal, a 
genetic marker of interest (e.g., a disease-associated marker), a common experience 
(e.g., traumatic stress/disaster, abuse victim, drug addiction), or a common learned 

30 ability (e.g., literacy), (e.g., compare juveniles with learning v. no learning). The 

subjects of the plurality may have the same gender, same age, within a 20, 15, 10, 5 or 
2 years. In one embodiment, each subject of the plurality is female, and the images 
were obtained at a similar phase of the menstrual cycle (e.g., the same quarter of the 
cycle). In one embodiment, each subject of the plurahty is addicted to a substance 

35 (e.g., a narcotic, caffeine). In one embodiment, each subject of the plurahty has a 
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5 abnonnal characteristic in a behavioral paradigm, e.g., a social reward paradigm, a 
CPT / probability paradigm, a physiological aversion / pain paradigm, a mental 
rotation paradigm, an emotional faces paradigm, and a monetary reward paradigm. 

[0116] In one embodiment, the virtual brain structure is smaller than normal, 

e.g., in one or both hemispheres. The brain structure can be the amygdala or other 
1 0 structure listed in Table 1 . 

[0117] The aligning can include locating one or more of the midpoints of 

decussations of the anterior and posterior commissures and the midsaggital plane. It 
can include rotation and/or a nondeformation transformation. 

[01 18] Deterniining positional information can include gray/white matter 

15 segmentation and/or evaluating signal intensity histograms, 

[0119] The methods can further include receiving information about the brain 

structure of an individual subject and comparing the received information to the 
information about the virtual brain structure. The method can fiarther include 
providing an estimate of risk for a behavioral trait, wherein the plurality of subjects 
20 each have a common behavioral trait. 

[0120] In another aspect, the invention features a datastructure that includes 

morphometric information about a brain structure (e.g., the amygdala or other brain 
structure). The morphometric information can be based on a statistical function 
dependent on a cohort of individuals with a common behavioral trait (e.g., an 

25 abnormal behavioral trait). For example, the morphometric information can include 
information about volume of the amygdala (e.g., a quantitative measure of volume) or 
information about the surface topology of the amygdala (e.g., information about the 
degree of undercutting or overcutting relative to a reference individual or a reference 
cohort or information about the surface contours, e.g., coordinates). In one 

30 embodiment, the information about surface topology of the amygdala describes at 
least a part of the right amygdala, or the right and left amygdala. In another 
embodiment, the morphometric information describes a degree of symmetry- 
asymmetry for a particular brain structure (e.g., the amygdala or other brain structure). 
A morphometric parameter described herein can be used as one parameter in a 
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5 classification method or any other method described herein, e.g., to evaluate 
genotypes and/or phenotypes and correlations between such factors. 

[0121] In another aspect, the invention features a method that includes: 

obtaining a group of subjects, e.g., human subjects; imaging the CNS of each subject 
while the respective subject is exposed to information (e.g., text, audio (e.g., music, 

_ • 

10 speech), video (e.g., advertising) etc); evaluating correlation between a characteristic 
of neural circuit activity of the subjects and alleles present at one or more genetic 
markers; and providing an evaluation of the information as a function between the 
characteristic and the firequency of an allele in a population. The method can include 
other features described herein. 

1 5 [01 22] In one exemplary method, subjects (e.g., human patients) are imaged 

using a pluraUty of procedures to produce tomographic maps. Generally, at least two 
(e.g., at least three, four, five, or six) different procedures are used. The plurality of 
procedures can include fimctional imaging during one or more paradigms, 
morphogenetic mapping of anatomical features, diffusion tensor analysis for white 

20 matter, radiological imaging, f-deoxy-glucose scanning, cerebral blood flow, and 
% cellular viabihty. 

[0123] Raw image data are translated into a multi-dimensional quantitative 

"systems biology map" that provides a complex representation of neuropsychiatric 
fimction. Because multiple procedures are used, the representation can span more 
25 than one cognitive center. Certain combinations of procedures can produce a nearly- 
continuous map that is a holistic measure of neuropsychiatric fimction. 

[0124] The systems biology map (SBM) can be displayed to a user as a matrix 

or may even be rendered on as a three-dimensional image of the brain. More 
typically, the SB map is stored in a database for computational analysis. Data can be 
30 analyzed using a models, e.g., to assess the reward-aversion circuit, e.g., using 
behavioral economics models. 

[0125] These SB maps have many apphcations, including, for example, 

evaluating a subject, diagnosing a subject, testing a therapeutic procedure or 
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5 therapeutic compound, monitoring disease progression, monitoring therapy, and so 
on. 

[0126] The fineness of the map may, for example, separate a general behavior 

perceived as a single disease into two or more distinguishable disorders. Further, the 
technique can be applied to non-human animals (e.g., primates and voles) and may be 
10 used in conjunction with administering a drug, evaluating gene expression, and so 
forth. 

[0127] The following are some exemplary features: a dataset that includes 

functional tomography for more than one paradigm; a dataset that includes parameters 
representing properties of more than one behavioral circuit; a multi-dimensional 
1 5 matrix that is a function of imaged nemral activity during a behavior and imaged 

anatomical features, (etc. for other combinations; a multi-dimensional matrix that is a 
function of at least three different images of the brain. 

[0128] An exemplary method of correlating a neuropsychiatric trait with a 

genetic locus may include: obtaining imaging information and genetic information 

20 from a population of individuals; generating a multi-dimensional systems biology 
(SB) map for each individual of the population; quantitatively sort the individuals 
based on their respective "maps", e.g., using an association rule algorithm, thereby 
identifying a subpopulation; comparing polymorphisms at least one genetic locus 
between individuals of the subpopulation to evaluate linkage between a 

25 polymorphism and members of the subpopulation 

[0129] For example, the comparing can include a genome scan to identify a 

genetic marker with a significant LOD score for the subpopulation. The method can 
also include comparing polymorphisms of individuals excluded from the 
subpopulation, e.g., to detect whether absence of an allele is determinative. Other 
30 genetic methods (e.g., families, linkage disequilibriiun, etc.) can be incorporated. 

[0130] After a genetic polymorphism is associated with a neuropsychiatric 

trait, a bottom up approach is used to evaluate individuals who have the 
polymorphism. The individuals can be evaluated at the extremes of function, and 
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5 imaged as described above to produce a SB map. Typically, the individuals that are 
evaluated are not members of the study that linked the polymorphism to the trait. 

[0131] This approach can have the following applications: provide 

confirmatory information, provide information for construction of a second model of 
neuropsychiatric function, and enable extrapolation of genetic information to a 
10 second population of individuals. 

[0132] The sorting can use criteria for at least two dimensions of the SB map. 

[0133] Many of the methods described herein can be embodied as software, 

e.g., a machine-executable instructions. The software can be stored on a machine- 
readable or accessible medium or as an article, e.g., a commodity. Such methods can 

15 also be implemented on a machine. Many steps within such methods can be executed, 
e.g., by interaction with a user or automatically. Methods can also be implemented 
across a network, e.g., an intranet or intemet. For example, the network can link a 
health care provider and a patient, a physician (e.g., a radiologist) and a patient, and 
different physicians (e.g., a radiologist and psychiatrist). Communications between 

20 members of the network can be secure, web-accessible, and can include hypertext, 
rotatable images, and other interactive and/or cartographic display techniques. 

[0134] The methods can be implemented using a system that includes a server 

that stores a database described herein, and a client that is in communication (e.g., 
digital comminication) with the server and can send queries and receive requested 

25 information to and firom the server. A server can include, for example, a processor 
and a memory for storing information about a plurality of subject. The processor can 
be configured to access the memory and retrieve information about one or more 
subjects and/or perform an analysis described herein. [0135] Some 
implementations of the invention enable providing a continuous function of disease 

30 risk. 

[0136] As used herein the term "circuit" refers to identifiable regions of the 

brain that are operational during a function such as a paradigm or other task. 
Typically such regions are distributed in space, but interact with one anottier. The 
brain is not modular but is a distributed system. 
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[0137] All cited patents, patent applications, and references are incorporated 

by reference in their entireties. In particular, U.S. published appUcations 2002- 
0042563 (09/822,585) and US applications 09/729,665, 60/573,138, and a provisional 
patent application, filed 2 August 2004, titled "MORPHOMETRIC ANALYSIS OF 
BRAIN STRUCTURES," bearing attorney docket number 00786-620P02, are 
incorporated by reference in their entireties. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0138] FIG. 1 is a schematic of exemplary interactions between the 

environment, genome, and epigenome. 

[0139] FIG. 2 is a schematic of an exemplary hierarchical organization that 

generates behavior. 

[0140] FIG. 3 is a schematic illustrating some levels of the components shown 

in FIG. 2. 

[0141] FIGs. 4a, 4b, and 4c are exemplary models for cognition. 

[0142] FIG. 5 is an exemplary model of the informational backbone for 

motivation (IBM) 20. 

[0143] FIG. 6 is an exemplary combined model 1 0 for motivation that depicts 

the interaction between iBM 20 and a behavioral mechanism 40 and a selection 
mechanism 30. 

[0144] FIG. 7 depicts the combined model for motivation 10 and interactions 

among its components. 

[0145] FIG. 8 depicts iBM components and mapping of five exemplary 

paradigms onto iBM components using a color code. 

[0146] FIG. 9 depicts aspects of an exemplary monetary reward paradigm and 

results obtained in particular experiments. 

[0147] FIG. 10a depicts an exemplary qualitative systems biology map. 

[0148] FIGs. 10b and 1 1 highhght some regions that feature in exemplary 

paradigms. 
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5 [0149] FIG. 12 depicts an exemplary qualitative systems biology map. 

[0150] FIG. 1 3 depicts brain regions which may feature in three putative 
endophenotypes for major depressive disorder. 

[0151] FIG. 14 depicts brain regions which may feature in some different 
disorders. 

10 [0152] FIG. 15 is a schematic of an exemplary phenotype -> genotype 
approach. 

[0153] FIG. 16 is a schematic of an exemplary genotype -> phenotype 
approach. 

[0154] FIG. 17 is a flow chart of an exemplary method for producing a QSBM 

4 

1 5 (quantitative systems biology map). 

[01 55] FIG. 1 8 is a flow chart of an exemplary method for generating a 

classification tree. 

[0156] FIGs. 19A, B, and C are exemplary methods for associating genotypic 

and phenotypic information. 

20 [0157] FIG. 20 depicts an exemplary set of matrices. 

[0158]. FIG. 21 is a schematic of an exemplary system 300, 

[0159] FIG. 22 depicts an exemplary computing unit. 

[0160] FIG. 23 is a schematic of an exemplary apparatus. 

[0161] FIG. 24 shows biiming across a spectrum. FIG. 25 are schematics of a 

25 binary tree. 

DETAILED DESCRIPTION 

[01621 The Informational Backbone of Motivation (IBM) 

[0163] One central features of the mind is the "informational backbone of 

motivation" or "IBM." The IBM is a large domain of the brain that processes 
30 information for motivation. See, e.g., FIG. 8. The IBM encompasses a number of 
circuits which participate in motivation. The IBM includes a number of component 
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mechanisms, including, e.g., mechanisms for representation and convergence, feature 
evaluation, probability assessment, outcome processing, valuation, reward/aversion 
processing, counterfactual comparisons, and memory. Paradigms can trigger one or 
more of these mechanisms, although not every paradigm evokes every circmt or 
structure in the iBM. [0 1 64] Other central features of the mind are depicted 

in FIG. 6. Exemplary component circuits include the reward/aversion circuit, 
working memory, centers of language and social behavior. . Still other components 
are mvolved in valuation, outcome processing, probabiUty assessment, feature 
evaluation, representation & convergence, reception, counterfactual comparisons, and 
other behaviors. It is possible to select or design paradigms that evoke one or more of 
these components and evaluate their function during the paradigm, e.g., as described 
for the exemplary paradigms provided herein. 

[0165] The reward/aversion circuit is part of the iBM. The reward/aversion 

circuit allows the organism to assign a value to objects in the environment so as to 
work for "rewards" and avoid "punishments" (aversive outcomes). This circuit can 
include an extended set of subcortical gray matter regions (nucleus accumbens (NAc), 
caudate, putamen, sublenticular extended amygdala (SLEA) , amygdala, 
hippocampus, hypothalamus, and thalamus) and domains of the paralimbic girdle 
(including the orbitofrontal cortex (GOb), insula, cingulate cortex, parahippcampus, 
and temporal pole) that receive dopaminergic projections from the ventral tegmental 
area and substantia nigra, here jointly refeired to as the ventral tegmentum: VT). 

[0166] Some additional exemplary regions of the brain that can be imaged and 

described in a systems biology map are provided in Table 1. 

Table 1: Exemplary Brain Regions 



1 . Transverse Cerebral Fissue and 
Third Ventricle 

2. Optic Chiasm 

3 . Fourth Ventricle 

4. Brainstem 

5. Lateral Ventricles 
6- Caudate 

7. Putamen 

8. Nucleus Accimibens 

9. Pallidum 



36. Occipital Pole 

37. Paracingulate Gyrus 

38. Precuneous Cortex 

39. Parahippocampal Gyrus, anterior 
division 

40. Parahippocampal Gyrus, posterior 
division 

41 . Parietal Operculum Cortex 

42. Postcentral Gyrus 

43 . Planimi Polare 
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10. 


Thalamus 


44. Precentral Gyrus 


11. 


Ventral Diencephalon 


45 . Planum Temporale 


12. 


Inferior Lateral Ventricles 


46. Subcallosal Cortex 


13. 


Amygdala 


47. Supracalcarine Cortex 


14. 


Hippocampus 


48. Supramarginal Gyrus, anterior 


15. 


Angular Gyms 


division 


16. 


Intracalcarine Cortex 


49. Supramarginal Gyrus, posterior 


17. 


Cingulate Gyrus, anterior division 


division 


18. 


Cingulate Gyrus, posterior 


50. Superior Parietal Lobule 


division 


5 1 . Superior Temporal Gyrus, anterior 


19. 


Cuneal Cortex 


division 


20. 


Central Opercular Cortex 


52. Superior Temporal Gyrus, 


21. 


Superior Frontal Gyrus 


posterior division 


22. 


Middle Frontal Gyrus 


53. Middle Temporal Gyrus, anterior 


23. 


Inferior Frontal Gyrus, pars 


division 


opercularis 


54. Middle Temporal Gyrus, posterior 


24. 


Inferior Frontal Gyrus, pars 


division 


triangularis 


55. ' Inferior Temporal Gyrus, anterior 


25. 


Frontal Medial Cortex 


division 


26. 


Frontal Operculum Cortex 


56. Inferior Temporal Gyrus, 


27. 


Frontal Orbital Cortex 


posterior division 


28. 


Frontal Pole 


57. Temporal Fusiform Cortex, 


29. 


Heschl's Gyrus (includes HI and 


anterior division 


H2) 




58. Temporal Fusiform Cortex, 


30. 


Insular Cortex 


posterior division 


31. 


Juxtapositional Lobule Cortex 


59- Middle Temporal Gyrus, 


(fonneriy Supplementary Motor Cortex) 


temporooccipital part 


32. 


Lingual Gyrus 


60. Inferior Temporal Gyrus, 


33. 


Occipital Fusiform Gyrus 


temporooccipital part 


34. 


Lateral Occipital Cortex; inferior 


6 1 . Temporal Occipital Fusiform 


division 


Cortex 


35. 


Lateral Occipital Cortex, superior 


62. Temporal Pole 


division 





5 



Each region of the brain can perform processes, and can be dedicated to a particular 
process. For example, the ventral prefrontal cortex is a communications center, e.g., 
in symbolic systems, e.g., language. Decision making is a behavioral output. 
Directed behavior can be the manifestation of intertwined processes of cognition, e.g., 
10 emotion and analysis. Various evaluation methods described herein can be used to 
identify deficits in one or more processes in one or more regions, thereby 
characterizing emotion and analysis,£0167] An Exemplarv Schema 
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5 [01 68] Referring now to FIG. 1 , behaviors can be explained by the interaction 

of three major factors, the genome, the epigenome, and the environment. The genome 
refers to the sequence content of nuclear genomic nucleic acid and mitochondrial 
genomic nucleic acid and other resident nucleic acid, such has viral nucleic acid. The 
epigenome refers to interaction of the genome with epigenetic factors that are 
1 0 transmissible, but variable modifications such as methylation, chromatin structure, 
long-range chromosomal effects (e.g., position effect variegation, transvection), and 
even RNAi (e.g., endogenous or exogenously added). The epigenome can function as 
a rheostat that reacts to create changes in biological function that are transmissible to 
a subsequent generation. 

15 [01 69] Referring now to FIG. 2, systems biology functions at the interface 

between behavior and the genome and epigenome. The genome and epigenome can 
directly affect cellular functions. At a higher level, groups of cells interact, e.g., as 
neural networks. At a still higher level distributed groups are formed which can 
control behavior, e.g., by reacting to the environment. Althou^ cells are critical 

20 component of the highest systems biology level, the impact of a single nucleotide in 
the genome or a single epigenetic factor can be only a very small jfraction of the 
complexity of the system. Thus, a single genetic or epigenetic change may be 
difficult detect in the noise of the system. 

[0170] Referring now also to FIG. 3, a variety of methods are available to 

25 obtain information about each level in the systems biology hierarchy. The 

information can include both structural and functional infomiation. For example, 
distributed neural groups can be evaluated by one or more of: magnetic resonance 
imaging (MRI) (also referred to as nuclear magnetic resonance or NMR) and other 
non-invasive techniques such as magnetic resonance spectroscopy (MRS), 
30 electroencephgraphy (EEG), magnetoencephalography (MEG), positron emission 
tomography (PET, including labeled hgand studies), optical imaging (OR), single 
photon emission computer tomography (SPECT), and functional computerized 
tomography (fCT). MRI methods include functional magnetic resonance imaging 
(fMRI), which provides information about function of neural groups. The map may 
35 include, for example, structural information such as moiphometric information about 
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5 anatomical features, diffusion tensor analysis for white matter, radiological imaging, 
f-deoxy-glucose scaiming, cerebral blood flow, and % cellular viability. 

[0171] Local circuits (e.g., neural groups) can be detected, e.g., by 

multicellular recording (e.g., during surgery of humans or by monitoring non-himians) 
or even by high resolution (or "fine") tomography, for example, by evaluating 50 |jM 

10 isotropic voxels at 7 T. Exemplary methods for evaluating cells include: evaluating 
local field potentials (LFPs), e.g., using implanted electrodes; evaluating ion or 
electrochemical potentials and flux (e.g., Ca^^ cascades, e.g., by voltometry); and 
evaluating gene and protein expression (e.g., using microarrays, antibodies, and mass 
spectroscopy). Methods for evaluating the genome and epigenome are described 

15 below. Many methods refer generally to genetic markers and genetic analysis. Such 
methods can also include evaluating epigenetic features associated with such markers. 

[01 72] Information firom different levels of the hierarchy can be combined. 

For example, mechanistic explanations can be derived by reductive hnkage of 

descriptions across (both up & down) scales. ' 

20 r01731 Strategies for Relating Phenotvpe and Genotvpe 

[0174] Three examples of general strategies for relating genetic markers to a 

phenotype defined by a systems biology map are described as follows. 

[0175] Referring to the example in FIG. 19A, the first strategy 110 includes 

first classifying 1 12 subject accordingly to their phenotype, e.g., using information 
25 jfrom systems biology maps (e.g., QSBMs). The classification process defines a 
plurality of phenotypic classes. Genetic markers are evaluated 114 for their 
association with (e.g., within) at least one of the classes. 

[0176] Referring to the example in FIG. 19B, the second strategy 120 includes 

first classifying 122 subject accordingly to their genotype, e.g., using genetic 
30 information. The classification process defines a pluraUty of genotypic classes. 
Then, phenotypes (e.g., information from systems biology maps) are evaluated to 
identify associations with at least one of the genotypic classes. 
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[0177] Referring to the example in FIG. 19C, the third strategy 130 includes 

concurrently classifying subjects by phenotype 1 12 and classifying them by genotype 
122. By exchanging information during the classification processes, a convergent 
solution can be obtained 1 34 that associates genotype and phenotype. For example, 
aconvergence of results can be forced between the neuroimaging and the genotypic 
data. This convergence relies on using outcomes firom the evalution of the 
neuroimaging data as a set of association rules to prune the partitions found with the 
data mining of the genotypic data. In parallel with this process, the outcome of the 
evaluation of the genotypic data is used to constrain the outcome from the 
neuroimaging data. 

[0178] Many aspects of the first strategy 110 - which involves classification 

by phenotype - have the further advantage of providing diagnostic and prognostic 
categories that are useful even without genetic information or validated genetic 
associations. Also, this first strategy does not necessarily require extensive family 
information, linkage disequiUbrium, foimder effects, or requirements on the mput 
population of subjects to provide meaningful statistics for finding genes that are 
associated with a particular phenotypic class. 

r01791 Producing an Exemplarv Systems Biologv Map 

[0180] Referring to the example in FIG. 8, a plurahty of paradigms can be 

used to generate a systems biology map that includes functional information about 
neural circuitry in a subject. 

[0181] FIG. 1 7 provides a flowchart for one exemplary method. A subject is 

evaluated using fMRI during a first paradigm 21 1 and during a second paradigm 212. 
Methods for evaluating subjects during paradigms are described, e.g., US 2002- 
0042563 and below. Raw acquisition data can be mapped 214 onto a standard 
anatomical model, e.g., the Talairach coordinates. In other implementations, other 
types information can be used instead of or in conjunction with the first and second 
paradigms. Such information includes: anatomical and morphological information 
about brain structure and function, clinical information (see, e.g., below). See 
discussion below. 
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[01 82] This information can be condensed 2 1 6 to produce a systems biology 

map, e.g., a qualitative or quantitative systems biology map. The abbreviation 
"QSBM" is used to represent quantitative systems biology maps. Although it is 
possible to use the raw data directly, typically the "raw" or "native" dataset acquired 
by instruments (e.g., an MRI machine) during a paradigm is very large. For example, 
a typically dataset from fMRI can include multiple 128x128 sections for 15 to 30 
different slices and for about 300 time points. If multiple runs of the paradigm are 
done, then the dataset is increased that many times. Parcellation and statistical 
analysis can further increase the dataset 10 to 15 fold. Thus, it is possible to have at 
least 20 Mb or even up to 1 terabyte (1 Tb) of data for a single subject. However, 
this information can be processed to generate a matrix (e.g., in the kilobyte to 1 Mbyte 
range) that has a reduced size relative to the native dataset, but retains the useful 
information. Thus, byte-for-byte, information can be condensed at least 10, 10^, 10^, 
10"*, 10^, or 10^ fold, and ranges therebetween. The ability to condense information 
into a meaningful and accessible format may be critical for the development and/or 
analysis of large databases of functional information about neural circuitry. 

[0 1 83] Extracting information for a QSBM typically involves discarding data. 

Although compression techniques can be used, e.g. to store the QSBM, typically the 
QSBM is in a form that is easily accessible, e.g., for computation. Because data is 
typically discarded, it is usually not possible to regenerate the native dataset from the 
QSBM. In one embodiment, the QSBM discards time resolution. For example, the 
QSBM can merely retain a list of activations for each region without reference to the 
temporal dimension, although the list may be ordered according to time of occurrence. 

[01 84] In one implementation the information is condensed into one or more 

matrices. FIG. 20 illustrates a set of matrices. Each matrix includes one dimension 
(illustrated vertically) that refers to different regions of the brain (e.g., regions 1, 2, 3, 
. . . n, wherein n refers to the n* region) and another dimension (horizontal) that refers 
to left and right hemispheres of the brain. A third dimension (e.g., going into the 
page) is used .to store a list of values for particular region-hemisphere. For fMRI data, 
for example, the Hst of values may refer to % change of each of the activation peaks 
detected during a paradigm. For example, if a region has three different activation 
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5 peaks that are due to the following changes in signal, 2.3%, -0.5%, and 1.2%, the 
sequence of values in the third dimension can be {2.3%, 1.2%, -05%}. 

[0185] Other values may also be used (e.g., instead of or in conjimction with 

% change), e.g., time to peak for each of the activation peaks detected during a 
paradigm, delay, FT, and slope. In another embodiment, information about each 
10 activation peak can further include information indicating location of the peak, e.g., 
where within the region the activation occurred and/or time, e.g., a reference to the 

r 

temporal dimension. 

[0186] An additional matrix can be used to store clinical (e.g., diagnostic) and 

demographic information such as age, gender, handedness, EEG, drug regimen (e.g., 

15 pharmacology), narcotic dependency, pedigree information, place of birth, place of 
residence, socio-economic status, race, language (e.g., ability to speak, understand a 
particular language and/or exposure to language, e.g., as an infant, child, adult), 
WAEB-R, DRM-rV diagnosis, and so forth. Still other types of useful information 
include quantitative medical assessments, e.g., blood pressure, pulse, body 

20 temperature, blood cell coimt, circadian rhythm, height, height, and other biometric 
values. 

[01 87] Another further matrix can be used to store genotypic information, 

although such information can also be stored separately. This additional matrix may 
only be two-dunensional. 

25 [0188] It is appreciated that a matrix can also be represented using other 

formats (e.g., an n-dimensional vector) or transformed into other representations (e.g., 
one or more tables in a relational database, a text string, and so forth). A set of 
matrices can also be represented as a single matrix which has an additional dimension 
relative to the most complex matrix in the set. 

30 [0189] FIG. 10 describes an exemplary qualitative systems biology map. The 

map describes different regions, here, the Gob, NAc, SLEA, Amygdala, and Vt. The 
map also indicates activity of the regions during the expectancy phase of cocaine and 
monetary reward paradigms, and the outcome phase of cocaine, monetary reward, 
beauty and pain paradigms. Using paradigms, component circuits or structures can be 
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5 evoked at one or more phases, e.g., during the expectancy or outcome phase. The 
expectancy and outcome phases refer to processes, hi some cases, these may occur 
concurrently, e.g., during overlapping time segments. 

[0190] Phenotvpic Classifications 

[0191] Referring now to the exemplary method 230 in FIG. 18, information 

10 about a plurality of subjects (e.g., human subjects) can be used to produce phenotypic 
classifications. The method 230 includes randomly selecting 232 a test set and 
training set from the plurality of subjects. For example, if a database includes 
information about 1000 subjects, 500 might be used for the training set and the other 
500 might be used for the test set. Variables for phenotypic information are then 
15 analyzed by evaluating 232 autocorrelations between the variables in the training set. 
The variables that are analyzed might include, e.g., all variables related to fimctional 
and structural information about neural circuitry. In other some embodiments, it may 
be usefiil to exclude the diagnostic and demographic information firom the 
autocorrelation analysis, although in other embodiments such information may be 
20 included. 

[0192] The results of the autocorrelations are then used to select variables to 

build 236 a classification tree. For example, the variable with the best autocorrelation 
score can be used to define a rule for the first node of the tree. The variable with the 
second best score may be used to define a rule for the second node of the tree. Details 
25 of tree building are provided below. Once the tree is complete, the classification is 
used to evaluate the test set. 

[0193] Objective scoring can also be used to evaluating the tree. For example, 

' the sizes of the clusters in the test set and the training set should be reasonably 
similar, e.g., within statistically acceptable values. In one embodiment, the tree 
30 should have a statistical significance, e.g., the tree structure is not attributed to chance 
alone. 

[0194] In some implementations, the classification method may achieve one 

or more of the following advantages: (i) classifications are obtained by completely 
objective criteria, (ii) structural and fimctional information can be easily integrated as 
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can other variables, e.g., information from different imaging techniques, different 
times and different subjects, (iii) classification is scalable and expandable (e.g., as the 
number of available subjects grows), (iv) information is condensed relative to raw 
acquisition data or data transformed onto an anatomical model. 

[0195] The use of autocorrelations enables one objective approach to selecting 

variables for tree classification. Variables that provide high autocorrelation scores are 
indicated as being highly informative. However, in some implementations, this 
objective approach is combined with a subjective approach, or if desired, in some 
implementations, a completely subjective approach can be used. In an exemplary 
subjective approach, regions of the brain that are known to connect and interact with 
regions that score high by objective criteria are also used for tree classification, e.g., 
independent of their own autocorrelation score. Regions that are knovm to be 
involved or be featured in a particular process may also be selected. 

[0196] Tree editing can include pruning branches, e.g., particularly branches 

that do not segregating individuals in an informative manner. For example, 
segregating a single individual firom a group of twenty does not aid the classification 
process. Similarly segregating in a upper node, five subjects from a group of five 
hundred may not inform the classification process. Pruning can be performed 
manually or automatically. In one example, associations rules are used to test the 
salience of possible correlations and to prune off non-informative nodes. In another 
example of automated pruning, branches with asymmetric distributions (e.g., < 10% 
into one branch) are removed by computer software. 

[01 97] The classification process can also be evaluated (e.g., by a user or 

software) to determine if it provides explanatory power. For example, a classification 
can be evaluated to determine if it bins known exophenotypes (e.g., clinical 
diagnoses) into subclasses. In another example, a classification is evaluated to 
determine if there is familiality, e.g., whether the classification identifies an 
endophenotype (see, e.g., below). In still another example, the classification 
continues until one or more particular constraints are satisfied. 

[0198] It is also possible to do a recursive process wherein tree branches are 

added and pruned diuing multiple recursive cycles until the tree structure satisfies 
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5 particular parameters, e.g., optimization parameters. For example, the tree can be 
modified until a cost value for growing the tree exceeds the informational value of the 
added complexity. 

[0199] In other embodiments, the training and test sets may be different sizes. 

Or in one embodiment, all available data is used to generate the classification tree. 

10 r02001 Endophenotvpes 

[0201] By evaluating information firom a plurality of subject it is possible to 

identify at least two types of phenotypic markers: endophenotypes and markers of 
. disease/disorder progression (MDP). 

[0202] An endophenotype typically includes the following properties; a) it 

15 ■ provides an internal marker of a probability Amotion for disease susceptibihty or 
resistance; (b) it is unchanged by illness progression; and (c) it has measurable 
heritability / familiality. See, e.g., Ahnasy and Blanquero (2001) Am, J. Med. Genet 
108:42. Thus, endophenotypes may be found (but not necessarily) in unaffected 
siblings andj)arents of a subject who is affected by a disorder. Similarly, the 
20 endophenotype can be present prior to onset of the disordCT. Thus, endophenotypes 
have high diagnostic value. An endophenotype may be defined by one or more 
variables, e.g., one or more variables present in a SBM (e.g., a QSBM) described 
herein. 

* [0203] In contrast, a marker of disease/disorder progression (MDP) is changed 

25 during the progression of a disorder. Such markers can be used to characterize the 

disorder, prescribe or monitor a treatment, and make other decisions (e.g., medical or 

financial decisions). 

[0204] A method for evaluating neuropsychiatric phenotypes can include a 

longitudinal component which is of great value in differentiating between 
30 endophenotypes and MDPs, Such longitudinal studies include analy^ng a subject at 
a first time and then analyzing the subject at a later time, e.g., at least one week, one, 
two, three, four, six, ten, or twelve months later. For example, the subject might be 
analyzed once a year over three to five years. In some embodiments, the subject is 
evaluated at approximately regular intervals. During these studies phenotypic 
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5 variables that remain unchanged, but which differ from normal (e.g., which are 

identified as useful for classification) are variables that can serve as endophenotypes. 
If the subject's outward clinical manifestations of a disorder are changing, other 
variables detected by evaluating neural circuit function may also change. Such 
variables can server as MDP. 

10 r02051 An Integrated Svstem 

[0206] Referring to FIG. 21, an exemplary integrated system 300 can be used 

to produce information for a database and generate information about neural circuit 
activity. For exiample, the system can include a network 305 that connects one or 
more imagers 350 (e.g., MRI machines) and one or more genotyping stations 340 with 

15 a database server 320. The imagers 350 can deliver raw or processed infonnation to 
the server 320 with infonnation that references an individual (e.g., using an 
anonymous index). The database server 320 also receives similarly referenced 
infonnation about the individual's genotype so that there is an association between the 
genotypic information and the phenotypic infonnation obtained by MRI. For 

20 example, a datastructure can be used that includes a first field with a pointer to the 
genotypic information of the individual and a second field with a pointer to the 
phenotypic information for the same individual. 

[0207] In one embodiment, the system 300 also includes a statistics engine 

which can evaluate the phenotypic information and/or genotypic infonnation, e.g., 
25 iising a method described herein. 

[0208] The methods and other features described herein can be implemented 

m 

in digital electronic circuitry, or in computer hardware, firmware, software, or in 
combinations thereof. Methods can be implemented using a computer program 
product tangibly embodied in a machine-readable storage device for execution by a 
30 programmable processor; and method actions can be performed by a programmable 
processor executing a program of instructions to perform fimctions of the invention by 
operating on input data and generating output. For example, methods can be 
implemented advantageously in one or more computer programs that are executable 
on a programmable system including at least one programmable processor coupled to 
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5 receive data and instructions from, and to transmit data and instructions to, a data 
storage system, at least one input device, and at least one output device. Each 
computer program can be implemented in a high-level procedural or object oriented 
programming language, or in assembly or machine language if desired; and in any 
case, the language can be a compiled or interpreted language. Suitable processors 

1 0 include, by v^ay of example, both general and special purpose microprocessors. A 

processor can receive instructions and data from a read-only memory and/or a random 
access memory. Generally, a computer will include one or more mass storage devices 
for storing data files; such devices include magnetic disks, such as internal hard disks 
and removable disks; magneto-optical disks; and optical disks. Storage devices 

15 suitable for tangibly embodying computer program instructions and data include all 
forms of non-volatile memory, including, by way of example, semiconductor memory 
devices, such as EPROM, EEPROM, and flash memory devices; magnetic disks such 
as, internal hard disks and removable disks; magneto-optical disks; and CD-ROM 
disks. Any of the foregoing can be supplemented by, or incorporated in, ASICs 

20 (application-specific integrated circuits). 

[0209] Data structures, trees, databases, and other information formats 

described herein can be stored in a machine accessible memory (e.g., volatile or non- 
volatile memory, within a CPU or external to a CPU) or on machine-readable medium 
(e.g., a hard disk, CD-ROM, and so forth. 

25 [02 10] An example of one such type of computer is depicted in FIG. 22, 

which shows a block diagram of a programmable processing system (system) 510 
suitable for implementing or performing the apparatus or methods of the invention. 
The system 510 includes a processor 520, a random access memory (RAM) 521, a 
program memory 522 (for example, a writable read-only memory (ROM) such as a 

30 flash ROM), a hard drive controller 523, and an input/output (I/O) controller 524 
coupled by a processor (CPU) bus 525. The system 510 can be preprogrammed, in 
ROM, for example, or it can be progranmied (and reprogrammed) by loading a 
program from another source (for example, from a floppy disk, a CD-ROM, or 
another computer). 
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5 [021 1] The hard drive controller 523 is coupled to a hard disk 530 suitable for 

storing executable computer programs, including programs embodying the present 
invention, and data including storage. The I/O controller 524 is coupled by means of 
an VO bus 526 to an I/O interface 527. The I/O interface 527 receives and transmits 
data in analog or digital form over conununication links such as a serial link, local 
1 0 area network, wireless hnk, and parallel link. 

[0212] , One non-limiting example of an execution enviroxmient includes 
computers running Linux Red Hat OS, Wmdows XP (Microsoft), Windows NT 4.0 
(Microsoft) or better or Solaris 2.6 or better (Sun Microsystems) operating systems. 
Browsers can be Microsoft Intemet Explorer version 4.0 or greater or Netsc^e 
15 Navigator or Communicator version 4.0 or greater. Computers for databases and 
^administration servers can include Windows NT 4.0 with a 400 MHz Pentium n 
(Intel) processor or equivalent using 256 MB memory and 9 GB SCSI drive. For 
example, a Solaris 2.6 Ultra 10 (400Mhz) with 256 MB memory and 9 GB SCSI drive 
can be used. Other environments can also be used. 

20 r02131 Diagnosis 

[0214] In one embodiment, a tree classification is produced based on 

information from a plurality of subjects. This tree can be used directly for diagnosing 
a subject (the "query subject"), particularly a subject that is not a member of the 
plurality of subjects that was used to produce the tree. Information for the query 
25 subject can be run through the tree. 

[0215] For example, if a native dataset for the query subject is received, the 

native dataset can be processed to produce a QSBM that has the same structure as the 
maps used for producing the tree. The query subject's QSBM is then compared to 
^ rules at each node of the tree to determine where the query subject falls on the tree, 
30 By proceeding down the tree to a terminal node, this process should indicate which 
bin or class the query subject belongs in. If the tree includes a probabiUstic or other 
statistical fimction that corresponds to the decision at each node, this process can also 
produce a probability or statistical significance for the diagnosis. For example, it is 
possible to display a value for each of the possible bins or classes that indicates the 
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5 probability that the query subject belongs in that bin or class. (The probabilities 
should sum to 1.0). However, it may not be necessary to explore all the branches of 
the tree. For example, only branches likely to be relevant might be tested. 

[0216] In another embodiment, a rules-based function is used to dejfine a class. 

Information about the query subject is then compared to one or more rules to produce 
10 an evaluation indicating whether the query subject belongs in the class. The result of 
the evaluation might again be a probability or other statistic. In this embodiment, it is 
not necessary to sequentially process a set of rules. 

[0217] Exemplary Applications 

[0218] There are numerous applications for the methods, data-structures, and 

15 systems described herein. In one example, the methods can be used to characterize 
(e.g., diagnosis) a neuropsychiatric disorder or a propensity or association with a 
disorder. In another example, the methods can be used for the discovery of a gene or 
epigenetic factor which contributes at least in part to a neuropsychiatric disorder. 
Such disorders include, e.g., schizophrenia, manic depression, bipolar disorder, 
20 addictions (e.g., substance abuse, gambling, etc.), obsessive-compulsive disorder, 
anxiety/paranoia, autism, schizo-affective disorder; delusional disorder, psychotic 
disorders not elsewhere specified; antisocial personality disorder, anorexia/bulimia 
nervosa; and so on. Similarly socially valued traits can also be evaluated, e.g., in 
individuals gifted with musical talent, charm, charisma, mathematical ability, 
25 persuasion, determination, creativity, and so forth. Once a gene or epigenetic factor 
is discovered it can be used a target for identifying, testing, or designing 
pharmacological interventions. 

[0219] Another exemplary application provides a database, which can be 

used, e.g., to diagnosis, evaluate, and process clinical or commercial information. 
30 The methods described herein can diagnose fiinctional brain disorder using multiple 
quantitative variables (e.g., by oversampUng information space. 

[0220] Still other applications include staging clinical diagnosis of 

neuropsychiatric in terms of functional impairment caused by non-pschiatric illness or 
to stage a psychiatric ilfaiess; detecting non-clinical variants that may appear as 
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5 clinical disorders; evaluating and planning treatment of psychiatric illness; monitoring 
and evaluating treatment efficacy; intervening in narcotics abuse; and monitoring 
narcotic consumption. [0221] Methods of Evaluating Genetic Material 

[0222] There are numerous methods for evaluating genetic material to provide 

genetic information. Genetic information can be obtained by evaluating a subject or a 

10 sample from a subject. The sample typically includes nucleated cells, e.g., somatic 
cells, or nucleic acid extracted from such cells (e.g., genomic DNA or cDNA or 
mRNA). In embodiments in which genomic DNA is used, virtually any biological 
sample (other than pure red blood cells) is suitable. For example, convenient tissue 
samples include whole blood, semen, saliva, tears, urine, fecal material, sweat, buccal, 

15 skin and hair. In embodiments in which cDNA or mRNA is used, the tissue sample 
usually includes cells in which the target nucleic acid is expressed. 

[0223] Nucleic acid samples can analyzed using biophysical techniques (e.g., 

hybridization, electrophoresis, and so forth), sequencing, enzyme-based techniques, 
and combinations-thereof. 

20 [0224] For example, hybridization to microarrays can also be used to detect 

polymorphisms, including SNPs. In one implementation, a set of different oligo- 
nucleotides, with the polymorphic nucleotide at varying positions with the 
oligonucleotides can be positioned on a nucleic acid array. The extent of 
hybridization as a function of position and hybridization to oUgonucleotides specific 

25 for the other allele can be used to determine whether a particular polymorphism is 
present. See, e.g., U.S. 6,066,454. 

[0225] In one implementation, hybridization probes can include one or more 

additional mismatches to destabilize duplex formation and sensitize the assay. The 
mismatch may be directly adjacent to the query position, or within 10, 1\ 5, 4, 3, or 2 
30 nucleotides of the query position. Hybridization probes can also be selected to have a 
particular Tm, e.g., between 45-60^C, 55-65°C, or 60-75**C. In a multiplex assay, 
Tm's can be selected to be within 5, 3, or 2*'C of each other, e.g., probes for a genetic 
marker can be selected with these criteria. 
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5 [0226] U.S. Pat. No. 5,837,832 describes a tiling method for array fabrication 

whereby probes are synthesized on a solid support. These arrays include a set of 
oligonucleotide probes such that, for each base in a specific reference sequence, the 
set includes a first probe (for example, a so-called "wild-type" or "WT" probe) that is 
exactly complementary to a section of the sequence of the chosen firagment including 

1 0 the base of interest in a first allele and at least one additional probes (called 

"substitution probe"), which are identical to the WT probe except that the base of 
interest has been replaced by one of a predetermined set of nucleotides (typically, 
one, two or three nucleotides), i.e., nucleotides other than the nucleotide in the first 
probe, for example a nucleotide complementary to a second allele. Probes may be 

15 synthesized to query each base in the sequence of the chosen firagment or a particular 
base known to be polymorphic. Target nucleic acid sequences which hybridize to a 
probe on the array which contain a substitution probe indicate the presence of a single 
nucleotide polymorphism. See also, e.g., U.S. 5,858,659; 5,861,242; 5,593,839 aad 
5,856,101 (describing, e.g., variously methods of using computers to design arrays 

20 and lithographic masks and methods of detecting insertions and deletions). 

[0227] The design and use of allele-specific probes for analyzing 

polymorphisms is described by e.g., Saiki et al.. Nature 324, 163-166 (1986); 
Dattagupta, EP 235,726, Saiki, WO 89/1 1548. Allele-specific probes can be designed 
that hybridize to a segment of target DNA firom one individual but do not hybridize to 

25 the corresponding segment from another individual due to the presence of different 
polymorphic forms in the respective segments firom the two individuals. 
Hybridization conditions should be sufficiently stringent that there is a significant 
difference in hybridization intensity between alleles, and preferably an essentially 
binary response, whereby a probe hybridizes to only one of the alleles. In one 

30 embodiment, probes are designed to hybridize to a segment of target DNA such that 
the polymorphic site aligns with a central position (e.g., in a 15-mer at the 7 position; 
in a 16-mer, at either the 8 or 9 position) of the probe. This design of probe achieves 
good discrimination in hybridization between different allelic forms. In one 
embodiment, the probes include a second mismatch which is non-complementary to 

35 both alleles of a biallehc pair. The second mismatch serves to destabilize the duplex, 
reduce Tm, and increase sensitivity. 
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[0228] Allele-specific probes are often used in pairs, one member of a pair 

showing a perfect match to a reference form of a target sequence and the other 
member showing a perfect match to a variant form. Several pairs of probes can then 
be immobihzed on the same support for simultaneous analysis of multiple 
polymorphisms within the same target sequence. 

[0229] Other hybridization based techniques include sequence specific primer 

binding (e.g., PGR or LCR); Southern analysis of DNA, e.g., genomic DNA; 
Northem analysis of RNA, e.g., mRNA; fluorescent probe based techniques (see, e.g., 
Beaudet et al. (2001) Genome Res. ll(4):600-8); and allele specific amphfication. 
Enzymatic techniques include restriction enzyme digestion; sequencing; and single 
base extension (SBE). These and other techniques are well known to those skilled m 
the art. 

[0230] Electrophoretic techniques include capillary electrophoresis and 

Single-Strand Conformation Polymorphism (SSCP) detection (see, e.g., Myers et al. 
(1985) ^ijft/re 313:495-8 and Ganguly (2002) i^MwMwto/. 19(4):334-42). Other 
biophysical mettiods include denaturing high pressure hquid chromatography 
(DHPLC). For example, different alleles can be identified based on the different 
sequence-dependent melting properties and electrophoretic migration of DNA in 
solution. Erhch, ed., PGR Technology, Principles and Applications for DNA 
Amplification, (W.H. Freeman and Go, New York, 1992), Chapter 7. Alleles of target 
sequences can also be differentiated using single-strand conformation polymorphism 
analysis, which identifies base differences by alteration in electrophoretic migration 
of single stranded PGR products, as described in Orita et al., Proc. Nat. Acad. Sci. 86, 
2766-2770 (1989). Amplified PGR products can be generated as described above, and 
heated or otherwise denatured, to form single stranded amplification products. Single- 
stranded nucleic acids may refold or form secondary structures which are partially 
dependent on the base sequence. The different electrophoretic mobilities of single- 
stranded amplification products can be related to base-sequence differences between 
alleles of target sequences. 

[023 1] In one embodiment, allele specific amplification technology that 

depends on selective PGR amplification may be used to obtain genetic information. 
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5 Oligonucleotides used as primers for specific amplification may carry the mutation of 
interest in the center of the molecule (so that amplification depends on differential 
hybridization) (Gibbs et al (1989) Nucleic Acids Res. 17:2437-2448) or at the 
extreme 3' end of one primer where, under appropriate conditions, mismatch can 
prevent, or reduce polymerase extension (Prossner (1993) Tibtech 1 1:238). See also, 

10 e.g., WO 93/22456. In one embodiment, the allele specific primer is used in 

conjimction with a second primer which hybridizes at a distal site. Amplification 
proceeds from the two primers, resulting in a detectable product which indicates the 
particular allelic form is present. A control is usually performed with a second pair of 
primers, one of which shows a single base mismatch at the polymorphic site and the 

15 other of which exhibits perfect complementarity to a distal site. 

[0232] ^ addition, it is possible to introduce a restriction site in the region of 

the mutation to create cleavage-based detection (Gasparini et al (1992) MoL Cell 
Probes 6:1). In another embodiment, amplification can be performed using Taq 
ligase for amplification (Barany (1991) Proc. Natl Acad. Sci USA 88:189). In such 
20 cases, ligation will occur only if there is a perfect match at the 3' end of the 5' 

sequence making it possible to detect the presence of a known mutation at a specific 
site by looking for the presence or absence of ampUfication. 

[0233] Enzymatic methods for detecting sequences include amplification 

based-methods such a the polymerase chain reaction (PGR; Saiki, et al. (1985) 

25 Science 230, 1350-1354) and ligase chain reaction (LCR; Wu. et al (1989) Genomics 
4, 560-569; Barringer et al (1990), Gene 1989, 1 17-122; R Barany. 1991, Proc, 
Natl Acad. iSci. C/iS4 1988, 189-193); transcription-based methods utilize 
synthesis by RNA polymerases to amplify nucleic acid (U.S. Pat. No. 6,066,457; U.S. 
Pat No. 6,132,997; U.S. Pat. No. 5,716,785; Sarkar et al.. Science (1989) 244:331- 

30 34; Stofler et al. Science (1988) 239:491); NASBA (U.S. Patent Nos. 5,130,238; 
5,409,818; and 5,554,517); rolling circle amplification (RCA; U.S. Patent Nos. 
5,854,033 and 6,143,495) and strand displacement amplification (SDA; U.S. Patent 
Nos. 5,455,166 and 5,624,825). Amplification methods can be used in combination 
with other techniques. 
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5 [0234] Other enzymatic techniques include sequencing using polymerases, 

e.g., DNA polymerases and variations thereof such as single base extension 
technology. See, e.g., U.S. 6,294,336; U.S. 6,013,431; and U.S. 5,952,174. For 
example, Chen et al., (PNAS 94:10756-61 (1997)), describes a locus-specific 
oligonucleotide primer labeled on the 5* terminus with 5-carboxyfluorescein (FAM). 

10 This labeled primer is designed so that the 3* end is immediately adjacent to the 

polymorphic site of interest. The labeled primer is hybridized to the locus, and single 
base extension of the labeled primer is performed with fluorescently-labeled 
dideoxyribonucleotides (ddNTPs) in dye-terminator sequencing fashion. An increase 
in fluorescence of the added ddNTP in response to excitation at the wavelength of the. 

15 labeled primer is used to infer the identity of the added nucleotide. 

[0235] Another method to identify SNPs is called single nucleotide primer 

extension (SnuPE) or minisequencing (Nikiforov et al.. Nucleic Acids Res., 22: 4167- 
75 (1994); Pastinen et al., Clin. Chem., 42: 1391-17 (1996); Landegren et al.. Genome 
Res., 8: 769-76 (1998); Kuppuswamy et al., Proc. Natl. Acad. Sci. U.S.A., 88: 1 143-7 

20 (1991)). This technique involves the hybridization of a primer immediately adjacent 
to the polymorphic locus, extension by a single dideoxynucleotide, and identification 
of the extended primer. All variable nucleotides can be identified with optimal 
discrimination using the same reaction conditions. (Pastinen et al.. Genome Res., 7:. 
606-14 (1997)). Related detection methods include luminous detection (Nyren et al., 

25 Anal. Biochem., 208: 171-5 (1993)), colorimetric ELISA (Nikiforov et al.. Nucleic 
Acids Res., 22: 4167-75 (1994)), gel-based fluorescent assays (Pastinen et al., Clin, 
Chem., 42: 1391-7 (1996)), homogeneous fluorescent detection (Chen et al.. Genet 
Anal., 14: 157-63 (1999)), flow cytometry-based assays (Cai et al.. Genomics, 66: 
135-43 (2000)), and high performance Uquid chromatography (HPLC) analysis 

30 (Hoogendoom et al.. Hum. Genet, 104: 89-93 (1999)). 

[0236] Mass spectroscopy (e.g., matrix assisted laser desorption ionization- 

time of flight (MALDI-TOF) mass spectroscopy) can be used to detect nucleic acid 
polymorphisms. In one embodiment, (e.g., the MassEXTEND™ assay, 
SEQUENOM, Inc.), selected nucleotide mixtures, missing at least one dNTP and 
35 including a single ddNTP is used to extend a primer that hybridizes near a 
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5 polymorphism. The nucleotide mixture is selected so that the extension products 
between the different polymorphisms at the site create the greatest difference in 
molecular size. The extension reaction is placed on a plate for mass spectroscopy 
analysis. See, e.g., Haff et al.. Genome Res., 7: 378-88 (1997); Griffin et al.. Trends 
Biotechnol., 18: 77-84 (2000); Sauer et al.. Nucleic Acids Res., 28: E13 (2000)). 

10 [0237] Fluorescence based detection can also be used to detect nucleic acid 

polymorphisms. For example, different temiinator ddNTPs can be labeled with 
different fluorescent dyes. A primer can be annealed near or immediately adjacent to 
a polymorphism, and the nucleotide at the polymorphic site can be detected by the 
type (e.g., "color^') of the fluorescent dye that is incorporated. 

15 [0238] It is also possible to directly sequence the nucleic acid for a particular 

genetic locus, e.g., by amplification and sequencing, or amplification, cloning and 
sequence. The direct analysis of the sequence can be accomplished, e.g., using either 
the dideoxy chain termination method or the Maxam— Gilbert method (see Sambrook 
et al., Molecular Cloning, A Laboratory Manual (2nd Ed., CSHP, New York 1989); 

20 Zyskind et al.. Recombinant DNA Laboratory Manual, (Acad. Press, 1988)). High 
throughput automated (e.g., capillary or microchip based) sequencing apparati can be 
used. In still other embodiments, the sequence of a protein of interest is analyzed to 
infer its genetic sequence. Methods of analyzing a protein sequence include protein 
sequencing, mass spectroscopy, sequence/epitope specific immunoglobulins, and 

25 protease digestion. 

[0239] Any combination of the above methods can also be used. For example, 

allele specific technology can be used in combination with microarrays. See, e.g., 
U.S. 6,287,778. 

[0240] Exemplary genetic markers (e.g., polymorphisms) can be found firom 

30 publicly available resources. Such resources include: the Whitehead Institute's 

integrated maps of the human genome (e.g., the WICGR map, Cambridge MA) which 
provide aligned chromosome m^s of genetic markers; other sequence tagged sites 
(STSs); radiation hybrid map data; CEPH yeast artificial chromosome (YAC) clones; 
the Genetic Annotation Initiative (web site: cgap.nci.nih.gov/GAI/; an NIH run site 
35 which contains information on candidate SNPs); dbSNP Polymorphism Repository 
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5 (world wide web site: ncbi.nlm.nih.gov/SNP/; a comprehensive NEH-run database 
containing information on SM*s and also haplotypes); HUGO Mutation Database 
Initiative (web site: arieLucs,imimelb.edu.au:80/.about.cotton/mdi.htm; a database 
with information about human mutations including SNPs); Human SNP Database 
(world wide web site: -genome.wi.mit.edu/SNP/human/index.html; managed by the 

1 0 Whitehead Institute for Biomedical Research Genome Institute, this site contains 
information about SNPs); SNPs in the Human-Genome SNP database (world wide 
web site: ibc.wustl.edu/SNP; providing access to SNPs that have been organized by 
chromosomes and cj^ogenetic location fromWashington University); HGBase (web 
site: hgbase.cgr.ke.se/; a sxmmiary of sequence variations in the human genome from 

15 the Karolinska Institute of Sweden); the SNP Consortiimi Database (web site: 

snp.cshl.org/db/snp/map; a collection of SNPs and related information resulting from 
a collaborative effort); GeneSNPs (world wide web site: genome.utah.edu/genesnps/; 
from the University of Utah and U.S. National Institute of Environmental Health). 
Many exemplary biallelic markers are also described in publications; see, e.g., U.S. 

20 Serial No. 60/206,615, U.S. Serial No. 60/216,745, WIPO Serial No. 

PCT/IBOO/00184, WIPO Serial No, PCT/IB98/01 193, PCX Publication No. WO 
99/54500, and WIPO Serial No. PCT/IBOO/00403, US 2002-0037508 and US 2002- 
0032319. 

[0241] The following are some examples of types of polymorphisms : A 

25 restriction fragment length polymorphism (RFLP) is a variation in DNA sequence that 
alters the length of a restriction fragment (Botstein et al.. Am. J Hum. Genet. 32, 314- 
331 (1980)). The restriction fragment length polymorphism may create or delete a 
restriction site, thus changing the length of the restriction fragment. RFLPs have been 
widely used in human and animal genetic analyses (see WO 90/13668; WO90/1 1369; 
30 Donis-Keller, Cell 51, 319-337 (1987); Lander et al.. Genetics 121, 85-99 (1989)). 
Other polymorphisms take the form of short tandem repeats (STRs) that include 
tandem di-, tri- and tetra-nucleotide repeated motifs. These tandem repeats are also 
referred to as variable number tandem repeat (VNTR) polymorphisms. VNTRs have 
been used in identity and paternity analysis (US 5,075,217; Armour et al., FEBS Lett. 
35 307, 1 13-115 (1992); Horn et al., WO 91/14003; Jefifreys, EP 370,719), and in a large 
number of genetic mapping studies. 
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5 [0242] Other polymorphisms take the form of single nucleotide variations 

between individuals of the same species. Such polymorphisms are far more jfrequent 
than RFLPs, STRs and VNTRs. Some single nucleotide polymorphisms (SNP) occur 
in protein-coding nucleic acid sequences (coding sequence SNP (cSNP)), in which 
case, one of the polymorphic forms may give rise to the expression of a defective or 

10 otherwise variant protein and, potentially, a genetic disease. Examples of genes in 

which polymorphisms within coding sequences give rise to genetic disease include P- 
globin (sickle cell anemia), apoE4 (Alzheimer's Disease), Factor V Leiden 
(thrombosis), and CFTR (cystic fibrosis). cSNPs can alter the codon sequence of the 
gene and therefore specify an alternative amino acid. Such changes are called 

IS "missense" when another amino acid is substituted, and "nonsense" when the 

alternative codon specifies a stop signal in protein translation. When the cSNP does 
not alter the amino acid specified the cSNP is called "silent". Other single nucleotide 
polymorphisms occur in noncoding regions. Some of these polymorphisms may also 
result in defective protein expression (e.g., as a result of defective spUcing). Other 

20 single nucleotide polymorphisms have no effect, e.g., no phenotypic effect. 

[0243] Pharmacologv and pharmacoeenomics 

[0244] It is also possible to use the methods described herein to evaluate 

phenotypes (e.g., by imaging) of a subject undergoing a treatment. Differences in 
phenotype can be detected by classification (e.g., classification trees). Then 

25 associations with a particular genotype can be detected. Other strategies (e.g., in FIG. 
19) can also be appUed, e.g., in combination with the data analysis methods and data 
structures described herein. Exemplary treatments include administering an agent 
(e.g., a medicament) and non-invasive treatments (e.g., hyponosis, psychotherapy, 
etc.). Homeopathic and traditional medicines as well as social behaviors can be 

30 similarly analyzed. 

[0245] In one embodiment, recursive partitioning is used in a study to do 

pharmacogenetics, e.g., using subjects undergoing a treatment (e.g., medication or a 
non-invasive therapy). Classification trees can be used to determine if subjects 
respond differently to a treatment. Or the classification can be done blind — e.g., 
35 evaluate treated subjects and controls to detect if significant classifications are 
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5 objectively made that discriminate between treated and untreated subjects (e.g., 
humans and non-humans.). 

r0246T Imaging 

[0247] An exemplary method for imaging a subject can include positioning 

10 subjects to be tested (e.g. persons who are under going a paradigm) and instructing the 
subjects to remain as still as possible, information about the brain is acquired. A 
measuring apparatus which non-invasively obtains information about the brain (e.g., 
structure and/or function) is used. In one embodiment, the subject to be tested is 
placed m a brain scanner, e.g., an MRI, flVIRI, MEG, fCT, OX, SPECT, or PET 
15 system. 

The imaged information can be acquired while the subject undergoes an experimental 
paradigm focused on one or more "motivation/emotion" processes. Altematively, 
signals can be acquired while the subject is exposed to certain stimuli (e.g. the subject 
views photographs of people or food or consumer products) or while the subject 
20 performs particular tasks (e.g. presses a bar to get a particular result). Altematively 
still, the subject can perform two or more of the above tasks while the CNS signals 
are obtained. 

The signals are statistically analyzed and localized to specific anatomical and 
functional brain regions. The details of the processes for statistically analyzing the 
25 CNS signals and localizing the signals to specific brain regions can vary. 

[0248] Referring now to the exemplary apparatus in FIG. 23, a noninvasive 

measurement apparatus and system for measuring indices of brain activity is 
described, e.g., as follows. In this particular example a magnetic resonance imaging 
(MRI) system 216 that may be programmed to non-invasively aid in the determination 
30 of indices of brain activity during motivational and emotional function in accordance 
with the present invention is shown. Its should be appreciated however that other 
techniques including but not limited to fMRI, PET, OI, SPECT, CT, fCT, MRS, MEG 
and EEG may also be used to non-invasively measure indices of brain activity during 
motivational and emotional function. 
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5 [0249] MRI system 215 includes a magnet 216 having gradient coils 216a and 

RF coils 216b disposed thereabout in a particular manner to provide a magnet system 
217. In response to control signals provided from a controller processor 218, a 
transmitter 219 provides a signal to the RF coil 216b through an RF power ampUfier 
220. A gradient amplifier 221 provides a current to the gradient coils 216a also in 
10 response to signals provided by the control processor 218. 

[0250] For generating a uniform, steady magnetic field required for MRI, the 

magnet system 217 may be provided having a resistance or superconducting coils and 
which are driven by a generator. The magnetic fields are generated in an examination 
or scamiing space or region 222 in which the object to be examined is disposed. For 
1 5 example, if the object is a person or patient to be examined, the person or portion of 
the person to be examined is disposed in the region 222. 

[025 1] The transmitter/amphfier combination 219, 220 drives the coil 216b. 

After activation of the transmitter coil 216b, spin resonance signals are generated in 
the object situated in the examination space 222, which signals are detected and are 

20 collected by a receiver 223. Depending upon the measuring technique to be executed, 
the same coil can be used as the transmitter coil and the receiver coil or use can be 
made of separate coils for transmission and reception. The detected resonance signals 
are sampled, digitized in a Digitzer/Aray proceser 224. Digitizer/Array processor 224 
converts the analog signals to a stream of digital bits which represent the measured 

25 data and provides the bit stream to the control processor 218. 

[0252] A display 226 coupled to the control processor 218 is provided for the 

display of the reconstructed image. The display 226 may be provided for example as a 
monitor, a terminal, such as a CRT or flat panel display. 

[0253] A user provides scan and display operation commands and parameters 

30 to the control processor 218 through a scan interface 228 and a display operation 

interface 230 each of which provide means for a user to interface with and control the 
operating parameters of the MRI system 215 in a manner well known to those of 
ordinary skill in the art. 
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5 [0254] The control processor 218 can be coupled to a signal processor 232 and 

a data store 236. The signal processor can be programmed according to a method 
described herein, e.g., to process raw image information. The processing can include 
localizing signals to a particular region of the brain. 

rQ255] Some Exemplary Brain Circuits 

1 0 [0256] Brain circuitry includes a prefrontal and sensory cortex. The prefrontal 

cortex includes medial prefrontal cortex and lateral prefrontal cortex. The region also 
includes the primary sensory and motor components. These components include the 
primary somatosensory cortex (SI), the secondary somatosensory cortex (S2) , the 
primary motor cortices (Ml), and secondary motor cortices (M2). Motor behavior 

15 involves regions such as Ml and M2, along with supplementary motor cortex (SMA). 
The frontal eye fields (102h) modulate motor aspects of eye control relating to 
directing the reception of visual signals from the environment to the brain. 

[0257] Brain circuitry also includes the thalamus region the dorsal striatum 

region and the lateral and medial temporal cortex regions. The medial temporal cortex 
20 region includes, for example, the hippocampus , the basolateral amygdala , and the 
entorhinal cortex. Also included as part of the brain circuitry are paralimbic regions 
which include, for example, the insula , the orbital cortex , the parahippocampus and 
the anterior cingulate. Current perspectives of reward circuitry also include the 
hypothalamus the ventral pallidum and a plurality of regions collectively designated. 

25 [0258] The regions collectively designated comprises the nucleus accumbiuns 

(NAc) the central amygdala the sublenticular extended amygdala of the basal 
forebrain SLEA/basal forebrain or SLEA/BF) the ventral tegmentum (ventral tier) and 
the ventral tegmentum (dorsal tier). 

[0259] The regions collectively represent a number of regions having i 

30 significant involvement in motivational and emotional processing. It should be 

appreciated that other components such as the basolateral amygdala are also important 
but not included in the regions designated by reference number. Other regions that are 
fiirther important to this type of processing include the hypothalamus, the 
orbitofrontal cortex , the insula and the anterior cingulate cortex . Further regions are 
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5 also important but listed separately such as the ventral pallidum , the thalamus , the 
dorsal striatum , the hippocampus , the medial prefrontal cortex , and the lateral 
prefrontal cortex . Not Usted in this figure but also involved in processing sensory 
infr)nnation for its emotional implications is the cerebellum. 

[0260] The functional contribution of each of these major regions are 

10 discussed below. It should be noted that what follows is a gross simplification and 
does not convey the complexity nor the diversity of the fimctions that these regions 
have been implicated with and may in the fiiture be connected to. Further note that 
there is currently a debate regarding the modular vs. non-modular fiinction of these 
brain regions, i.e., can a specific function be attributed to each region in isolation. 
1 5 Accordingly what is listed below is information which provides one of ordinary skill 
in the art with the understanding that this function may be mediated by the connection 
of this region with many other regions (i.e., the function mediated by a distributed set 
of regions, of which the identified region is a fundamental component). 

[026 1 ] As a brain region the NAc has previously been impUcated in the 

20 processing of rewarding/addicting stimuli, and is thought to have a number of 

functions with regard to probability assessments and reward evaluation. It has also has 
been implicated in the moment by moment modulation of behavior (e.g., initiation of 
behavior). Signals measured from the NAc are shown and described below in 
conjunction with FIGS. 3A-3D. 

25 [0262] The SLEA/BF has been imphcated in reward evaluation, based on its 

likely role in brain stimulation reward effects. It is thought to be important for 
estimating the intensity of a reward value. It and other sections of the basal forebrain 
appear to be important for the processing of emotional stimuU in general, and it has 
been implicated in drug addiction. 

30 [0263] like the NAc, the amygdala has been implicated in.both processing of 

emotional information along with processing of pain and analgesia information. The 
amygdala has been implicated in both the orienting to and the memory of 
motivationally salient stimuli across the entire spectrum from aversion to reward. It 
may be important for the processing of signals with social salience in real time. In this 

35 context it is often referred to with regard to fear. A number of its anatomical 
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5 connections to primary sensory cortices, suggest that it is important for the 
modulation of attention to motivationally salient stimuli. 

[0264] With respect to the VT/PAG, doparminnergic projections are present 

from the VT to the SLEA, tfee orbitofrontal cortex, the amygdala, and the NAc. 
Lideed dopaminergic projections go to most subcortical and prefrontal sites. In FIG. 

10 3, the fundamental importance of the VT/PAG projection is focussed on the NAc, 

central Amygdala, and SLEA/B , though it also projects to other regions. The VT has 
been implicated in reward prediction processes, motor functions and a number of 
learning processes around motivational events in general. The PAG has also been 
implicated as a modulator of pain stimuli, for example, and may therefore be a region 

IS that signals early information on rewarding or aversive stimuh. 

[0265] The GOb component of the prefrontal cortex has been imphcated in a 

number of cognitive, memory, and planning functions aroxmd emotional stimuli or 
regarding rewarding or aversive outcomes in animal and human studies. This section 
of the prefrontal cortex has also been implicated in modulating pain. It has afferent 
20 and efferent connections with a number of subcortical structures. The GOb is 

involved in a number of different reward processes including those of expectancy 
determination and reward valuation. Patients with lesions in this region tend to have 
impulse control problems. 

[0266] The hypothalamus is involved in the monitoring and maintenance of 

25 homeostatic systems. It also has been both implicated in the evaluation of the 
relevance for rewarding and aversive stimuli in order to maintain homeostatic 
equilibrium. The hypothalamus is highly important for meeting the objectives which 
optimize fitness over time and meet the requirements necessary for survival. 

[0267] The cingulate cortex has been interpreted to be involved in attention 

30 and plaiming, the processing of pain impleasantness, the processing of reward events 
and emotions in general, and the evaluation of emotional conflict. The cingulate 
cortex is an extensive region of brain cortex and appears to have emotional and 
cognitive subdivisions, to name a few. 



-55- 



wo 2005/020788 



PCT/US2004/025146 



5 [0268] The insula has been implicated in number of functions including the 

processing of emotional stimuU, the processing of somatosensory functions , and the 
processing of visceral function. 

[0269] The thalamxis is composed of a number of sub-nuclei which have been 

implicated in a diverse range or functions. Fundamental among these functions 
1 0 appears to be that of being an informational relay of sensory and other information 
between the extemal and intemal environment. It has also been directly implicated in 
both rewarding and aversive processes, and damage to.the structure may resuh in 
dysfunction such as chronic pain. 

[0270] The hippocampus has been extensively implicated in functions for 

1 5 encoding and retrieval of infonnation. Lesions to this structure lead to severe 
"impairment in the abihty to form new memories. Motivated behavior is heavily 
dependent on such memories: for instance, how a particular behavior in the past led to 
obtaining a goal object which would reduce a particular deficit state such as thirst. 

[0271] The ventral palladium is one of the primary output sources of the NAc 

20 and has a number of projection sites including the dorsomedial nucleus of the 

thalam\xs. Via this connection, it is one of the major relays between the NAc and the 
rest of the brain, in particular prefrontal cortical regions. It has been strongly 
implicated in reward functions and is a site thought to be important for the 
development of addiction. 

25 [0272] The medial prefrontal cortex of the brain has been strongly implicated 

in reward functions and has been found to be one of the few brain sites into which 
cocaine self-administration can be initiated in animals. 

[0273] In response to reward and aversion situations, certain regions of the 

brain circuitry play a role in processing reward/aversive information to plan 
30 behavioral responses as discussed above. These regions are designated 

reward/aversion regions of the brain. The activation of such reward/aversion regions 
can be observed during positive and negative reinforcement using neuroimaging 
technology. These reward/aversion regions produce specific functional contributions 
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5 to motivated behavior. For example, contributions made by regions such as the 
include assessment of probability. 

r02741 Morphometric information 

[0275] Morphometric information about brain structures provides useful 

indicators of resistance/susceptibility to a disorder or abnormal behavior, therapeutic 

10 efficacy, and the presence and staging of active illness (e.g., an active disorder or 

abnormal behavior). Morphometric information includes information that describes a 
spatial and/or structural property of a brain structure or relevant part thereof. An 
example of a brain structure is the amygdala. In the case of cocaine addiction, the 
right amygdala and certain subnuclei (e.g., as mentioned below) may be particularly 

15 relevant. Other examples of brain structures are provided in Table 1 . 

[0276] Morphometric information can be in the form of a moiphometric 

parameter, e.g., a quantitative or qualitative parameter. A quantitative measure of 
volume is one form of a morphometric parameter. Coordinates or equations 
describing a surface of a brain structure are. another example. Morphometric 

20 information can be absolute, e.g., relative to a particular coordinate-firame, or can be 
relative. For example, a linear distance measure of the extent of over- or under- 
cutting of a brain structure surface of a subject relative to a reference brain structure is 
a useful form of relative information. To illustrate, we have observed in one set of 
cocaine addicted individuals that the volume of the right amygdala is decreased about 

25 23% and that relative to an iso-surface based on probability 0.5, addicted individuals 
have an undercutting in the anterior extent of about 4.5 mm. 

[0277] The use of morphometric information, e.g., as described herein, will 

aid in the diagnosis of behavioral disorders and neuropsychiatric disorders as well as 
in the discovery of dmgs and other therapies for treating such disorders. 

30 [0278] In one implementation, a virtual reference structure is created, e.g., 

representing a statistical function for a brain structure among a cohort of individuals, 
e.g., individuals with a common characteristic. For example, the cohort can be a 
cohort of nonnal controls, a cohort of disorder affected individuals, e.g., substance 
affected individuals, or bipolar disorder-affected individuals. Using images taken of 
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S the brain or regions thereof, brain structures can be segmented as individual 
structures, following standardized anatomic definitions. Caviness, V. S., Jr., 
Kennedy, D. N., Richebne, C, Rademacher, J., and Filipek, P. A. (1996). The human 
brain age 7-11 years: a voliunetric analysis based on magnetic resonance images. 
Cereb Cortex 6, 726-736; Makris, N., Meyer, J. W., Bates, J. F., Yeterian, E. H., 

10 Kennedy, D. N., and Caviness, V. S. (1999). MRI-Based topographic parceliation of 
human cerebral white matter and nuclei n. Rationale and apphcations with 
systematics of cerebral connectivity. Neuroimage 9, 18-45; Makris, N,, Worth, A. J., 
Sorensen, A. G., Papadimitriou, G. M., Wu, O., Reese, T. G., Wedeen, V. J., Davis, T. 
L., Stakes, J. W., Caviness, V. S., et al. (1997). Morphometry of in vivo human white 

1 5 matter association pathways with diffusion-weighted magnetic resonance imaging. 
Ann Neurol 42, 951-962; Seidman, L. J., Faraone, S. V., Goldstein, J. M., Goodman, 
J. M., Kremen, W. S., Toomey, R., Tourville, J., Kennedy, D., Makris, N., Caviness, 
V. S., and Tsuang, M. T. (1999). Thalamic and amygdala-hippocampal volume 
reductions in first-degree relatives of patients with schizophrenia: an MRI-based 

20 morphometric analysis. Biol Psychiatry 46, 941-954; Seidman, L. J., Faraone, S. V., 
Goldstein, J. M., Kremen, W. S., Horton, N. J-, Makris, N., Toomey, R., Kennedy, D., 
Caviness, V. S., and Tsuang, M. T. (2002). Left hippocampal volume as a 
vulnerability indicator for schizophrenia: a magnetic resonance imaging 
morphometric study of nonpsychotic first-degree relatives. Arch Gen Psychiatry 59, 

25 839-849. Segmentation can be performed manually, semi-automatically, or 
automatically. 

[0279] Images can be registered (aligned), e.g., to a reference brain that was 

separate from the cohort. A probabihty surface (or "isoform surface" or "iso- 
surface") for a particular structure for each cohort can be calculated, e.g., on a voxel- 
30 by-voxel basis with the aligned data. Iso-surfaces for a pre-selected probability value 
(e.g., probability 0.5) are created for each cohort separately. Three-dimensional 
visuaUzation of these siufaces can be used to look for systematic differences in the 
topology of the brain structure between cohorts, or to evaluate a brain structure of a 
subject in comparison to the cohort. 
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5 [0280] The following non-limiting example illustrates some aspects of the 

methods described herein in one particular implementation 

[0281] Example (Part 1) 

[0282] Cocaine and nicotine are two of the most acutely reinforcing drugs in 

humans and in animals; they are also profoundly addicting, and have a strong co- 

1 0 morbidity with depression. Twenty-five percent of the U.S . Population suffers firom 
nicotine dependence, and smoking leads to about 500,000 deaths per year. Major 
depression is the most coramon psychiatric disorder in the U.S. today, and the number 
one cause of mortaUty in the world. It is frequently co-morbid in individuals that 
cease nicotine or cocaine self-administration. Baseline anhedonia or dysthymia is 

1 S also hypothesized to be a causal factor in the development of nicotine dependence and 
is observed in individuals between episodes of cocaine self-administrations. Long- 
term use of psychostimulants and the ensuing dependence has pronounced effects on 
the circuitry of reward-aversion. Recent neuroimaging and post-mortem stereology 
have also documented functional and morphometric changes in the brain circuitry of 

20 reward-aversion with mood disorders (Manji JSK, Drevets WC, Chamey DS, (2001) 
The cellular neurobiology of depression. Nat Med. 7:541-7.; Ongur D, Drevets WC, 
Price JL.( 1998) Glial reduction in the subgenual prefirontal cortex in mood disorders. 
Proc Natl Acad Sci USA. 95:13290-5.). 

[0283] The neural circuitry that mediates the rewarding (i.e., hedonic) effects 

25 of psychostimulants (Koob et al., 1998), or the rewarding and aversive effects of other 
stimuU (Wise et al., 1978; Wise, et al., 1992; Koob, 1992; Stein & Fuller, 1992; 
Kometsky & Esposito, 1981), can be readily studied by fMRI BOLD to obtain 
quantitative measures of changes in brain activity. These circuits include the: nucleus 
accumbens (NAc), sub-lenticular extended amygdala (SLEA) of the basal forebrain, 
30 amygdala, ventral tegmentum (VT), and orbital gyrus (GOb), along with other 
parahmbic regions such as the anterior cingulate, insula, parahippocampus, and 
temporal pole. Together, these brain regions (referred to as the reward-aversion 
circuitry) appear to be fundamental to the assessment of motivationally salient 
informational features for organizing behavior. 
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5 [0284] In humans, these bram regions have been shown to process expectancy 

and valuation information and the sequential effects of expectancy on subsequent 
outcomes. The differential valuation of rewarding vs. aversive outcomes utilizes the 
same brain circuitry, and unique signal profiles in a subset of these regions have been 
mapped for rewarding vs. aversive stimuli. We can characterize the relative 
10 contributions made by each of these subregions to discrete components of reward- 
aversion function in different individuals using paradigms (e.g., paradigms described 
in Breiter et al,, 1996; Cohen et al., 1996; Seidman et al., 1998; Breiter & Rosen, 
1999; Aharon et al., 2001; Becerra et al., 2001; Breitei: et al., 2001). Thus, we can 
sample, for example: 

IS (1) stimulus input and representation, 

(2) feature extraction necessary for assessing motivational intent in 

others, 

(3) probability functions necessary for expectancy determination , 

(4) expectancy vs. outcome functions, 
20 (5) valuation functions, and 

(6) positive vs. aversive outcomes. 

We can develop a systems biology map, e.g., by using information from at least two 
of these Auctions. The map can describe how a stimulus that is rewarding or aversive 
is processed. 

25 [028S] In this exemplary case, the system assessed determines reward- 

aversion function and acts as an informational backbone for motivation. The ability to 
produce such system biology maps in individuals further gives us a precise 
mechanism by which to characterize malfunctions in this circuitry that quantitatively 
characterize functional brain disorders such as stimulant addiction and depression. 

30 [0286] Circuitry-based events responsible for behavior and intracellular 

signaling events, at very different spatiotemporal scales of brain function, are 
interlinked and that processes at the distal ends of this spatiotemporal continuimi can 
serve as markers, e.g., for genetic analysis. 
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5 [0287] Analysis of brain structure and/or function produces a set of 

quantitative indices (e.g., a systems biology map) which can be associated with 
genetic information from the subject. Typically such genetic information includes 
markers on a plurality of different non-homologous chromosomes. When sufficient 
number of individuals are analyzed, statistics can be used to evaluate the relationship 
10 between genotype and phenotype. Linkage from a set of quantitative indices, such as 
the multitude of quantitative measures in a systems biology map, to the quantitative 
measures of molecular genetics can pinpoint the genes that contribute to susceptibility 
and/or resistance to frmctional brain disorders such as addiction and depression. 

[02881 Example CPart 2V. Detailed Methods 

15 [0289] (a) Subject recruitment, screening, and scheduling 

[0290] For Phi, a total of 500 subjects (pl\is 8% of this number as potential 

replacements) will be recruited for scanning over one year (months 6 - months 18 of 
the project, and then rescanned during months 19 - 30). For Ph2, a total of 800 
subjects per year (plus 8% of this mmiber as potential replacements) will be recmited 

20 over 4 years (total = 3200 + replacements). All phases of this project will be 

conducted according to the U.S. Food and Drug Administration guidelines and the 
Declaration of Helsinki. To protect all sensitive data, we will obtain a Certificate of 
Confidentiahty from NIH. Written informed consent will be obtained from all patients 
before protocol-specified procedures are carried out. Subjects will be drawn from an 

25 outpatient sample, and will be recruited through general media, as well as physician 
referrals. 

[0291] Inclusion Criteria'. The following conditions must be met for patient 

eligibility: 

(1) Written informed consent. 
30 (2) Men and women aged between 20-65 years, as sib pairs who are 

concordant for the criteria below, discordant, or in nuclear families with these 
diagnoses). 

(3) Nicotine dependent subjects: 

a. Smokers who have smoked>10 cigarettes/day for more than 2 years 
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b. Meet DSM-IVR criteria for Nicotine Dependence, as determined 
with the Fagerstrom Nicotine Tolerance Questionnaire (FTQ) (Fagerstrom, 
1978) 

c. Saliva cotinine levels > 14 ug/L, and end-expiratory carbon 
monoxide levels > 8ppm (subjects with alcohol dependence will be 
excluded). 

(4) Cocaine dependent subjects: 

a. DSM-IVR diagnosis of cocaine dependence (who are actively 
using cocaine at the time of entry and are not seeking or participating in 
treatment for addiction) without other Axis I psychiatric illnesses or past 

< 

experience of violent behavior while abusing cocaine or opiates. Exceptions 
will be made for dependence on caffeine or for mild (less than 3 standard 
drinks/week) consumption of alcohol. 

b. Validation of subjects self-reported drug use will be performed 
using hair specimens assayed for levels of commonly abused drugs. Our 
previously published data indicate that self-reported substance use in our 
non-treatment seeking research subjects is generally valid (Elman et al., 
1999). 

(5) Subjects with recurrent major depressive disorder: 

a. DSM-IVR criteria for lifetime diagnoses of Unipolar Depressive 
Disorders (major depression, dysthymia, and minor depression), according 
to the Structured CUnical Literview for DSM-IV - Axis I Disorders/Patient 
Edition (SCID-I/P) (First et al., 1995) and Inventory for Depressive 
Symptomatology (Rush et al., 1986, 1996; GulUon et al., 1998). The 
DSM-rVR SCID approach will be used by clinicians, and complemented 
with a SSAGA-II performed by trained research assistants. 

(6) Subjects with (3) and/or (4) and/or (5) (see Dierker et al, 2002) 

(7) Healthy controls without (3) or (4) or (5) 

[0292] Exclusion Criteria: la brief, subjects with any of the following will be 

excluded: pregnancy; suicidality or homocidality; serious medical illness including 
HIV + status; severe respiratory compromise; current use of nicotine-containing 
products in subjects without nicotine dependence; history of seizure disorder; 
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5 delirium, dementia, or mental disorders due to general medical conditions; substance 
abuse not specified above; schizophrenia; schizo-affective disorder; delusional 
disorder; bipolar disorder; psychotic disorders not elsewhere specified; antisocial 
personality disorder, unless comorbid with cocaine dependence; current 
anorexia/bulimia nervosa; clinical laboratory evidence of 
1 0 hypothyroidism/hyperthyroidism. 

[0293] (b) Experiments 

[0294] The six experimental paradigms listed below will be run with all 

subjects in Phi and Ph2. For each subject, these six paradigms will be run in the order 
in which they are listed. The time needed for each paradigm wiU be: 

15 #1 social reward paradigm, 22 minutes, 

#2 CPT / probability paradigm, 4 1/2 minutes, 
#3 physiological aversion / pain paradigm, 4 1/2 minutes, 
#4 mental rotation paradigm, 8 min and 54 seconds, 
#5 emotional faces paradigm, 11 1/2 minutes, 

20 #6monetary reward paradigm, 24 minutes. 

[0295] Between each paradigm, 2 minutes are scheduled for the overall 

imaging session to allow the reading of the next set of instructions. The total time for 
functional imaging will thus be 75 1/2 minutes, plus approximately 10 minutes for 5 x 
2 minute pauses between the scanning of each paradigm. This thus leaves 35 minutes 
25 for structural scanning described in the imaging section. Of these 6 paradigms, 4 of 
them have a traditional block design (#2 - #5), while 2 of them (#1 & #6) have a 
single trial-like design. These paradigms have been chosen because they robustly 
activate reward-aversion circuitry, to produce a systems biology map. 

[0296] (1) Social reward paradigm (Aharon et aL, 2001) 

30 [0297] Social stimuli will consist of two sets of 40 non-famous human faces 

[digitized at 600 dpi in 8-bit grayscale, spatially down sampled, and cropped to fit in 
an oval "window" sized 310-350 pixels wide by 470 pixels high using Photoshop 4.0 
software (Adobe Systems)]. Each set will consist of 20 male and 20 female faces. 



1 
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5 Subjects will be told tbat they will be exposed to a series of pictures that if not 

interfered with, will change every eight seconds. However, if they want a picture to 
disappear faster, they can alternate pressing the "z" and "x" keys, whereas if they want 
a picture to stay longer on the screen, they can alternate pressing the "n" and "m" 
keys. The dependent measures of interest will be the amount of work in units of key 
10 press that subjects exert in response to the different categories of stimuh, and their 
resulting viewing durations. 

[0298] Each pair of key presses will be set to increase or decrease the total 

viewing time according to the following formula: NewTotalTime = OldTotalTime + 
(ExtremeTime-OldTotalTime) / where ExtremeTime was 0 seconds for keypresses 

15 reducing the viewing time, ExtremeTime will be 14 seconds for keypresses increasing 
the viewing time, and K was a scaling constant set to 40. If the elapsed time for the 
picture surpassed the total time determined by keypressing, the picture was removed 
and the next trial began. A "slider" was displayed left of each picture indicating total 
viewing time at any moment, and changing with every keypress. Subjects will be 

20 informed that the task will last 40 minutes, and that this length is independent of their 
behavior during the task, as is their overall payment for participating in the 
experiment. 

[0299] (2) CPT with differential probability conditions (Breiter & Rosen, 

1999; Seidman et al., 1998) 

25 [0300] The set of experimental conditions in this study are designed to parse 

out differences in vigilant attention during a serial processing continuous performance 
task [CPT-AX(del)], involving a simple probabilistic relationship between a cue and 
delayed target, versus a dual processing continuous performance task [CPT-AX(int)], 
with a complex probability relationship between a cue and delayed target. The 

30 conditional probability of a subsequent target, given the incidence of a cue, will be the 
same between tasks since the CPT-AX(del) and CPT-AX(int) tasks have the same 
total number of cue-target pairs, and the same total incidence of true cues plus false 
cues. The tasks will be different in that the determination of cue-target pairs is more 
effortful for the CPT-AX(del) task, due to divided processing and interference 
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5 suppression needs. The effortftil determination of cue-target pairs will impair 
probability computation and lead to dnninished task performance. 

[0301 ] The two paradigm conditions will involve computer presentation of an 

auditory letter string, with each letter spoken at a rate of 1 per second. These 
paradigms will have an A-B-A-C-A-C-A-B-A desijgn where the A condition will be a 

10 simple CPT (referred to as the "QA" sequence), and the B condition will be an 

effortful CPT with three letters between cue and target pairs. The B and C conditions 
will involve either serial processing (CPT-AX(del)) or divided/dual processing (CPT- 
AX(int)). The CPT-AX(del) is characterized by a lack of false cues or targets between 
each cue ("q") and target ("a") pair, or by any interdigitated cue-target pairs (i.e., 

15 "q'L"q'L"a"_"a"). thus allowing simple probabilistic assessment of cue to target 
pairing with serial association of stimulus and response. The CPT-AX(int) has false 
cues and/or targets between pairs of cues and targets, and has cue-target pairs 
interdigitate together so that commingled pairs were possible, thus preventing simple 
counting or rehearsal procedures (i.e., forcing subjects to maintain two or more 

20 counts), and increasing the effort needed for probabilistic assessment of cue to target 
pairing. Each B and C epoch will last 90 seconds, while the baseline A epochs will 
last 60 seconds. There will be a target to distracter ratio of .13 for both A and B 
conditions, and the nxmiber of cue-target pairs will be the same. Subjects will respond 
with a magnet compatible button press, so that reaction time and accuracy could be 

25 recorded. The order for performing the CPT-AX(del) and the CPT-AX(int) will be 
counterbalanced across subjects. 

/ 0302] (3) Physiological Aversion (Thermal Pain) (Becerra et aL, 2001 ) 

[0303] Subjects will be informed in detail about the nature of the experiment, 

and the temporal sequence of procedures, including rating methods. These rating will 

30 involve rating on a scale from 0 (no pain) to 10 (maximum pain) their perception of 
the pain they experienced, after the functional run. Thermal stimuli will be delivered 
using a modified [Becerra et al., 1999] Peltier based thennode (Medoc, Haifa, Israel). 
One scan will be performed during which a base temperature of 35 ""C (30 s) 
(condition A), a warm stimulus of 41 °C (25 s) (condition B), and a target temperature 

35 of 46°C (condition C) will be interleaved. The thermode will be set to change the 
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5 temperature at a rate of 4 °C/s. Thus, it will take 2 s to reach 41 °C from the baseline 
and 2 s to return to baseline, while for the 35-46°C contrast, the delay will be 4 
seconds. The delays were not part of the baseline (30 s) or stimulus (25 s) times. The 
three stimuh will be interleaved in a block design: A-B-A-C-A-C-A-B-A. 



They will thus consist of three-dimensional perspective drawings of 10 cubes 
arranged in chiral patterns and viewed from a variety of rotation angles. Two task 
variants will be used. In a control condition, subjects will be shown a pair of figures, 
half of which are identical, and half of which are miiror-reversed shapes. Each of the 

15 10 possible angled-shapes (0-1 80** in 20** increments) will appear in each type of 
pair. The stimulus ordering will consist of a set of blocks, so that each of the stimuli 
appear once before any stimulus appears twice, and each appears twice before any 
appears three times, and so forth. Within each of these blocks, the stimuli will appear 
in random order except that the same stimulus will not appear twice within three 

20 successive trials. Moreover, half of the pairs within each block will include identical 
figures and half include mirror-reversed figures. No more than three consecutive trials 
can have the same response. 

/ [0306] The second version of the task (a rotation variant) will be identical to 

the first except that the members of each pair will be presented at different 

25 orientations. The left member will always be presented so that the major axis is 

vertical. The right member will be presented at nine possible angles (20 - 1 80° in 20° 
increments) from vertical. Three sets of these rotation trials will be used (and 4 sets of 
control trials), which will include rotations aroimd different major axes. One set will 
include rotations aroimd the x-axis, another around the y-axis, and another aroimd the 

30 z-axis. These stimuli will be presented in separate sets. Within each set, the stimulus 
trials will be ordered so that each orientation appears once before it appears again, 
once with identical stimuli and once with mirror-imaged stimuli, within each balanced 
subgroup of 18 trials. The same orientation will not appear twice within three 
consecutive trials. 



f0304J 



(4) Mental rotation (Cohen et aL, 1996) 



10 [0305] 



The figures will be the original Shepard and Metzler (1971) objects. 
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5 [0307] A third !'resting" or fixation condition will be interleaved between the 

"control'* and ^'rotation" tasks. Subjects will be asked to look at each pair, and to 
decide whether the figures are identical or are nairror-images and to indicate their 
choice by pressing one of two buttons. In the control condition, subjects will be asked 
to simply respond as quickly and accurately as possible. In the rotation condition, they 
10 will be told to visualize the right-hand stimulus rotating until it is aligned with the 
left-hand stimulus, and then to decide whether the two shapes are identical or are 
mirror reversed. 

[0308] (5) Emotional faces (Breiter et oL, 1996) 

[0309] Faces used in these experiments will be firom Ekman and Friesen 

15 (1976). They will have been standardized by (a) digitization, (b) scaling of extents, (c) 
normalization of contrast across all expressions for each of the individuals utilized 
(N=8), and across all individuals in the cohort., and (d) fitting with an oval mask to 
minimize the observation of hair. 

[03 1 0] The experiment will employ an A-B-A-C-A-D- A-B-A-D-A-C-A-B-A 

20 design with equal length epochs of tachistoscopic-hke presentations of the faces as. In 
A, subjects will see 1 80 presentations of 8 faces in random order; neutral expressions 
(200msec) will be followed by a fixation point (300msec). In B, C, and D, subjects 
will see faces with one emotion presented 180 times per epoch with the same timing 
parameters as in A. These facial expressions will be: happy (B), angry (C), and 
25 fearfiil (D). The order in which these blocks of facial expression are presented will be 
counterbalanced by emotion, and by epoch order within run. This will be a covert 
paradigm design with passive viewing of tachistoscopic-like face presentations, and 
use of the same 8 individuals for each expression presented in random order per 
epoch. 

30 [0311] (6) Monetary expectancy, gains, and losses (Breiter et oL 2001) 

[03 12] In this experiment we seek to map the hemodynamic changes that 

anticipate and accompany monetary losses and gains under varying conditions of 
controlled expectation and counterfactual comparison. The display wiU consist of 
either a fixation point or one of 2 disks ("spinners"). Each spinner will be divided into 
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5 2 sectors. Both spinners will offer the same outcomes, a gain of $+10 or a loss of $-8, 
but the likelihood of the gain will be high (0.66) on the "good" spiimer and low (0.33) 
on the "bad" spinner. The relative areas of the spinner assigned to the two outcomes 
represent the likelihoods. Thus, on the good spiimer, 66% of the area is colored green 
and labeled $+10, and the remaining 33% of the area is colored red and labeled $-8; 
10 on the bad spinner, the colors and labels are reversed. Providing larger gains than 
losses will be implemented to compensate for the tendency of subjects to assign 
greater weight to a loss than to a gain of equal magnitude. 

[0313] Before the game begins, subjects will be shown each spinner 3 times 

so as to leam its composition. Bach trial will consist of (1) an "expectancy phase," 

1 5 when a spinner is presented and an arrow spins around it, and (2) an "outcome phase" 
when the arrow lands on one sector and the corresponding amoimt is added to or 
subtracted from the subject's winnings. During the expectancy phase, the image of one 
of the 2 spinners will projected for 10 sec, and the subject will score their emotional 
response to the displayed spinner (or fixation point) using a potentiometer. During 

20 the outcome phase, the arrow will land on one of the sectors and flicker for 9.5 

seconds, indicating how much they won or lost. During this time, subjects will score 
their emotional response to the observed outcome. After 9.5 seconds, a 0.5 second 
mask will appear. On fixation-point trials, an asterisk will appear in the center of the 
display for 19.5 sec, followed by the 0.5-sec notask. The pseudo-random trial sequrace 

25 will be fiilly counter-balanced to the first order so that trials of a given type (spinner + 
outcome) are both preceded and followed by the same number of 4 spixmer/outcome 
combinations and 2 times by fixation-point trials. Subjects will observe 24 trials of the 
+$10 outcome, 24 trials of the -$8 outcome, and 16 trials of spinner baseline. A 
"dummy" trial will be inserted at the beginning and end of each run for 

30 counterbalancing, allowing 18 trials per run for 4 runs. Runs will be separated by 2 
min rest periods. The same trial sequence will be used for all subjects, generating 
winnings of $48, to which will be added the $50 endowment. 

[0314] (c) Imaging (3T and 7T) 

[0315] Five hundred subjects in Phi and potentially 3200 subjects in Ph2, plus 

35 replacement subjects, will be scanned on a 3.0 T AUegra System (Siemens) using a 
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S quadrature Siemens head coil. The Siemens system performs a whole head shimming 
procedure before scanning begins, which incorporates a full array of second order 
shims to optimize BO homogeneity, and thus reduce susceptibility/resistance in 
targeted reward-aversion regions of interest. Imaging for all experiments will begin 
with a 3-plane scout scan (conventional FLASH sequence with isotropic voxels of 2.8 

10 mm). The axial and coronal scouts will be used for placement and prescription of a 
3D MPRAGE anatomic scan, which will be used for anatomic locahzation of 
functional activation, and quantitative volumetric measurements. Prescription of 
experimental slices will follow this sequence with 30 shces parallel to the AC-PC line 
and covering the NAc, amygdala, SLEA, hypothalamus, VT and GOb, along with 

IS most of the lateral prefrontal cortex, and components of the parietal-occipital junction. 
BOLD imaging will then be performed using a gradient echo T2* weighted sequence 
(TRyTE=2000/29 ms,; FOV = 20 cm; in-plane resolution 3.125x3.125 mm, slice 
thickness = 3 mm; 30 contiguous axial slices). 

[0316] For Phi, 100 subjects will be scanned on a 7.0T ultra high field system 

20 developed for functional brain studies. If the results of comparison between the 3T 
and 7T systems are favorable to the 7T system, then the 3200 subjects of Ph2 will be 
scanned on it. The 7T system consists of a whole body magnet (Magnex Scientific) 
with a custom made resistive shim set (through 3rd order) and custom head gradient 
set. The study will obtain a 3D MPRAGE anatomic scan, and then a conventional T2 
25 scan at the same sUce locations as the functional prescription. Functional imaging wiU 
consist of a high resolution (1 .5mm x 1 .5mm x 3mm) single shot gradient echo 
sequence, covering less brain volume than the 3T scaiming protocol, but including the 
NAc, amygdala, SLEA, hypothalamus, VT and GOb, along components of the lateral 
prefrontal cortex and the parietal-occipital junction. 

30 [031 7] (d) Data analysis of neuroimaging data 

[0318] . Anatomic segmentation/parcellation for volumetrics and activation 
localization 

[0319] The anatomic scans of all subjects will undergo segmentation and 

parcellation. Segmentation methodology based on intensity contour and differential 
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5 intensity contour concepts can be used (see, e.g., Kennedy et al., 1989; Caviness et 
al., 1996; Filipek et aL, 1994; Rademacher et al., 1992). The cortical parcellation 
technique is based upon the concept of limiting sulci and planes and takes advantage 
of the observed relationships between cortical surface features and the location of 
functional cortical areas. A critical advantage of this method is that the definitions are 
10 unambiguously definable in a standardized fashion from the information visible in 
high resolution MRI. 

[0320] To perform this process with 500+ subjects for Phi and 3200+ subjects 

for Ph2, we can use an automated, fully 3D procedure for whole-brain segmentation. 
The technique uses a set of manually labeled brains as a training set in order to 

15 compute prior probabilities and class statistics, and applies a Bayesian classification 
rule. Specifically, we compute the maximum a posteriori (MAP) estimate of the 
segmentation W given an input image / and prior information firom the training set. 
Formally this can be expressed as maximizing p(W\I), the probability distribution of 
the segmentation given the observed image intensities. The prior probability of a 

20 given segmentation is initially encoded assuming that the classification at each voxel 
is independent of all other voxels. This constraint is then relaxed, and the image is 
iteratively resegmented using an anisotropic Markov-random field to model the image 
segmentation, resulting in a fimal segmentation that is more spatially uniform as well 
as more acciirate than the initial one. 

25 [032 1 ] Manually parcellated surfaces will also be used as a training set that 

can be employed to constmct classifiers in an analogous manner to the sub-cortical 
segmentation procedure. This process will depend on two properties of the cortical 
surface. The first is mean curvature of the surface, a differential measure of the 
surface folding computed from the trace of the Hessian matrix of the height function 

30 of the surface over its tangent plane at each point. The second is the average 

convexity of the surface (Fischl et al., 1999), a measure that is more sensitive to the 
presence of primary folds than to secondary or tertiary folds. The initial labeling is 
performed by assumiag the classification is spatially independent, so that the 
probability of the neuroanatomical label at each point in the cortex is independent 

35 from all other cortical locations. This is of course not the case, as the probability of 
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each label is related to the labels of the neighborhood in which it Ues. In order to 
capture the spatial regularity of the labehng, we model the surface labeling using an 
anisotropic Markov random field. The anisotoropy comes from the observation that 
labels are much more likely to change as one moves across the cortex in the direction 
of maximum curvature (i.e. the first principal direction) than in the direction of 
minimum curvature (i.e. the second principal direction). This information is encoded 
in the form of Gibbs priors on the probability of a given labeling. The most probable 
labeling is then iteratively recomputed using the independent spatial labeling as input, 
usmg the Iterated Conditional Modes (ICM) algorithm (Besag, 1974). 

[0322] To further investigate the surface-based structure of each subjects' 

brain, we will further use a set of automated tools for the construction of 
geometrically accurate and topologically correct models of the cortical sheet. These 
include accurate segmentation of gray matter and white matter (Dale et al., 1999), 
inflation and flattening of the surface models for visualization and analysis purposes 
(Fischl et al., 1999, 2000, 2002; Sereno et al., 1995), and automatic correction of 
topological defects (Fischl et al., 2002). The explicit construction of both the 
gray/white and pial surface boundaries allows the accurate measurement of the 
thickness of the cortical sheet (Fischl et al., 2000). The thickness of cortex is a 
potentially important diagnostic measure for a variety of neurodegenerative and 
psychiatric disorders, many of which are associated with progressive, regionally 
specific atrophy of the gray matter (for instance, consider alterations in prefrontal and 
temporal cortex volume observed with cocaine dependence; Franklin et al., 2002). 

I 

[0323] FMRI data preparation 

[0324] For this project, four of the experimental paradigms have a traditional 

block design, and 2 have a single trial-like design. For the block design experiments, 
data preparation will generally follow the procedures specified in Aharon et al. 
(2001), while for the single trial-hke designs, it will generally follow procedures 
specified in Breiter et al. (2001). In Phi, data preparation and assessment of main 
effects between groups (below) will involve analysis using FSL / FS Fast. As an 
example of the process planned for data preparation, data preparation will involve 
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motion collection, intensity normalization, signal detrending, and spatial filtering. For 
example, after motion correction, time series data will be inspected to ensure that no 
data set evidences residual motion in the foma of cortical rim or ventricular artifacts > 
1 voxel. Functional data will then be intensity scaled on a voxel-by-voxel basis to a 
standard of 1000, so that all mean baseline raw magnetic resonance signals are equal. 
These data will then be detrended to remove any hnear drift over the course of the 
scan. Spatial filtering will be performed using a Hanning filter with 1.5 voxel radius 
(this approximates a 0.7 voxel gaussian filter). Lastly, the mean signal intensity for 
each voxel over all runs will be removed on a time point by time point basis. For the 
single trial-like experiments, data will fiirther be selectively averaged and normalized 
relative to the 4 time points of data preceding the trial (see Breiter et aL, 2001). 

[0325] The data analytic procedures used in this project will be based on two 

assumptions, (a) that the behavior of the hemodynamic control system is 
approximately linear (i.e., it obeys the superposition axiom) under the experimental 
conditions tested and in the brain regions targeted by these paradigms, and (b) that 
deviations firom hemodynamic stationarity Mali be correctable by means of the 
normalization procedures employed. If the hemodynamic control system obeys 
superposition and stationarity, then the counterbalancing procedure used in each 
paradigm ensures that any carryover of hemodynamic responses fi-om antecedent 
experimental conditions will be constant across conditions. 

[0326] Two separate approaches will be applied to the evaluation of saUent 

changes related to experimental condition. The first will be based on the evaluation of 
individual data on its native anatomy. The second will be based on the evaluation of 
aggregate effects on averaged data that are then evaluated within all individuals in the 
cohort, using a standardized anatomical space. 

[0327] Individual data on native anatomy: Individual analyses will be pursued 

since aggregate analyses may produce Type I errors m the case of (1) opponent 
responses to different experimental conditions, which would tend to cancel as a result 
of averaging, or (2) responses confined to a small proportion of trial types or confined 
to a putative phenotype, which may be diluted by averaging. For single trial-like 
experiments, data obtained at all time points for each experimental condition will be 
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Statistically evaluated by correlation with a model impulse function (Boynton et al., 
1996; Dale & Buckner, 1997). To eliminate cross-trial hemodynamic overlap, 
statistical maps will be derived from correlation between the y function and a 
difference signal between each experimental condition and the paradigm baseline. For 
block design experiments, a y function will be convoluted with the experimental time 
course, and used in a correlation analysis. For botii single trial-like experiments and 
block design experiments, the outcome of correlation analysis will be assessed for 
foci of signal change using a cluster-growing algorithm (for example: Bush et al., 
1996). Clusters selected for further analysis will be required to either meet a corrected 
statistical threshold, or have signal intensity changes from baseline > 0.05%. For the 
corrected statistical threshold (in order to maintain an overall a < 0.05), the cluster- 
growing algorithm will localize activations that meet a corrected p value threshold of 
p < 0.00075 (0.05/67) for the number of segmentation and parcellation regions 
searched. Regions of interest (ROIs)will be dehneated by the voxels with p < 0.05 in a 
7mm radius of the voxel with the minimum p value (the "max vox"). These ROIs will 
then be used to sample the % signal intensity change per condition from the 
experimental baseline. 

[0328] Identified ROIs will be localized via superposition of the segmentation 

and parcellation contours produced as described above. The % signal intensity change 
from baseline for each of the experimental conditions in tiie task will then be 
quantitated and organized in a matrix based on the anatomic segmentation and 
parcellation units vs. hemispheric laterality. If a focus of signal change is observed in 
an anatomic segmentation/parcellation unit, it will be noted in the matrix for that 
anatomic region as an itemset of % signal changes from baseline of each experimental 
condition. Each experiment will have an independent matrix, as will the volumetric, 
and clinical information. 

[0329] Aggregated individuals in Talairach space: Analysis of individuals 

may miss low-level signal changes observed in common across the cohort, hence an 
analysis on aggregate data vnW also be utilized. The outcomes of this analysis that are 
not found by the analysis of individuals on native anatomy (above), will supplement 
the results found above. Indeed, these results will constitute a second matrix for each 
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5 experimental paradigm, thus producing a total of twelve matrices with functional data, 
along with two more from volumetric and clinical data. 

[0330] Analysis of aggregated individuals will identify foci of signal change 

across the aggregate, and apply ROIs of these foci to individuals to sample the % 
signal intensity change per condition from the experimental baseline. It must be noted 
10 that such analyses from the aggregate may produce Type I errors in the case of 

opponent responses to different experimental conditions, which would tend to cancel 
as a result of averaging, or responses confined to a small proportion of trial types or 
confined to a putative phenotype, which may be diluted by averaging. 

[0331] To allow averaging of data across subjects (for Phi, 250+ subjects and 

15 then 500+ subjects; for Ph2, 1600+ subjects, and then 3200+ subjects), each 
individual's set of functional data and structural data will be transformed into 
Talairach space (Breiter et al., 1996a, c; Talairach & Toumoux, 1988), and resUced in 
the coronal orientation with isotropic voxel dimensions (x,y, z = 3.125 mm for 3T, 
and =1.5 mm for 7T). Optimized fit between functional data and structural scans will 

20 then be obtained via translation of exterior contours. These Talairach transformed 
functional and structural scans will then be averaged. The same procedures for ROI 
identification used with the individual data on native anatomy will then be used to 
identify a set of activation clusters on the averaged data to then be used as ROIs to 
sample from each Talairach transformed individual data set the % signal intensity 

25 change per condition from the experimental baseline. In the averaged data, only 
activations will be selected that meet a corrected p value threshold of p < 0.00075 
(0.05/67) for the number of segmentation and parcellation regions searched. 
Activation ROIs from the averaged data will be anatomically localized in each 
individual on the Talairach transformed individual structural data, using superimposed 

30 segmentation and parcellation contours that have also been morphed into the 

Talairach domain. Again, as with the data produced by analysis of individuals, the 
data produced fix>m ROIs determined on averaged data, will be listed in matrices. 
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[0332] Classification tree analysis (phenotyping) 

[0333] To partition the nexiroimaging data into the fewest number of sets for 

the quantitative indices measured, that will be predictive of any future data set 
obtained, an algorithm based on classification tree analysis will be used. These 
analyses will be performed on the functional and quantitative volumetric data 
organized in matrices for each individual. In general, these techniques split data sets 
presented to them into sub-classes, and keep track of how it was done via a decision- 
tree structure. This decision tree structure can then be used to classify novel data. 
There are a number of classification techniques, all of which aim to select the class 
with the highest estimated conditional probability without computing the whole 
probabiUty distribution. These techniques basically differ m their employment of 
different biases in their first steps. Regression trees are basically like classification 
trees, with the difference that they handle continuous data, and given that the 
functional and structural data will be continuous, they will be utilized. For these 
algorithms, most of the effort goes into determining the optimal ordering of the 
variables in the decision tree, as well as the level at which to cease the decision 
process (i.e., there comes a point when all members of a branch should be in the same 
classification). These algorithms are also typically "non-parametric" in that no 
predictive model has to be hypothesized for the fitting. 

[0334] The software that will be used for this process will be the CART 

(classification and regression trees) system initially designed by Steinberg and Colla 
(1995) and distributed by Salford Systems, CA (note: this system is also incoiporated 
into many statistical packages such as S-plus). CART can lead to "over-fitting" of the 
data, in that it finds too many classes. Overfitting leads to the identification of too 
many itemsets (e.g., interesting patterns in the data), which can be a serious issue in 
domains with many multi-valued parameters, where the search space is large. To 
protect against this outcome, association rules (Srikant & Agrawal, 1995) can be used 
to test the salience of all the possible correlations between subclasses in the data, and 
then prune off non-informative decisions. A standard ^proach for making optunal 
class predictions using association rules is the "Large Bayes Classifier" of Meretakds 
and colleagues. In general, these techniques are computation intensive, necessitating 
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5 the use of a commercial cluster box system or supercomputer. A recmrsive- 
partitioning technique (see, e,g.,Zhang and Bonney, 2000) can also be used. 

[0335] During Phi , the first 60 subjects per diagnostic category scanned on 

the 3T magnet (total of 250 subjects) will be used as a training set, and the subsequent 
subjects scanned will be used as a test set to assess the initial classification schema. 
10 This process will then be repeated using the fiill complement of subjects scanned on 
the 3T as a training set, and the 100 subjects scannjed on the 7T as a test set. A greater 
specification of the identified classes foimd firom the larger training set would indicate 
that the initial cohort had not produced saturation of the identified endophenotypes. 
The training set size can be enlarged as the project progresses. 

15 [0336] Assessment of main effects between groups 

[0337] To evaluate the efiBcacy of our phenotyping methods, statistical 

assessment of effects between groups and correlations between fimctional and 
structural measures will be performed. These analyses should produce results 
embedded in the output of the classification tree analysis. Estimates of the central 

20 tendency (location) and dispersion (scale) of the data distribution of the diagnostic 
groups will use conventional least squares statistics or a robust statistics module, e.g., 
a Tukey bisquare estimator (Hoaglin et al., 1983; Breiter et al., 2001). Robust 
statistics are less subject than conventional parametric statistics to the influence of 
outUers and provide more eflScient estimates of location and scale of contaminated 

25 normal distributions. Although robust methods are more efficient when dealing with 
contaminated distributions, they are less efficient than parametric statistics when 
dealing with near-normal distributions. 

[0338] Main effects between experimental conditions, and experimental 

conditions in each diagnostic group can be assessed using multiple regression to carry 
30 out a random effects analysis of variance (ANOVA). Experimental condition will be 
defined as a categorical (noncontinuous) variable, thus avoiding any assumptions 
concemiag the form of the time courses. For the data determined on individual native 
anatomy, the ANOVA results will need to meet a more stringent a level than the 
conventional 0.05 value by correction for the number of clusters tested in each 
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5 individual. For the data determined from averaged data, the ANOVA results will need 
to meet a more stringent a level than the conventional 0.05 value by correction for the 
number of clusters found on the averaged data. For both data measured from 
individual native data, and data from Talairach transformed data, in cases that meet 
the criterion a level, pair-wise contrasts between specific experimental conditions will 
1 0 then be performed. 

[0339] As a last analysis, an autocorrelation analysis will be performed among 

the ftinctional and structural measures within diagnostic groups, and between 
diagnostic groups, using a Pearson product-moment correlation coefficient. The 
correction for performing multiple autocorrelations will be 0.05 adjusted by the 
1 5 nimiber of calculations performed. 

[0340] FMRI power analysis 

[0341] To estimate our statistical power to detect a difference between 

experimental conditions that would segregate potential endophenotjpes for cocaine 
dependent, nicotine dependent, or mood disordered subjects vs. controls, we first 

20 determined the expected effect-size by reanalyzing prior experimental data. We 
reanalyzed data from a set of experimental stimuli that produced similar signal 
magnitude changes in reward regions to those produced by the 6 experimental 
paradigms described above, and that had a similar number of time points to those 6 
paradigms in their shortened format planned for this project: We also selected data for 

25 these calculations from experiments that involved subjects with cocaine addiction, 
and subjects who were healthy controls. In one case we reanlyzed a prior cocaine 
infiision study (Breiter et al., 1997), while in another we reanalyzed a morphine 
infrision study in healthy controls (Breiter et al., 2000). For each subject in the 
cocaine infiision study, signal from all voxels in the bilateral NAc (defined 

30 anatomically on Talairach-transformed images) was normalized, averaged, and 
linearly detrended. The resulting time series of 136 time points had an average 
standard deviation, across 20 infusions in 13 subjects, equal to 0.84% of the grand 
average signal level. The difference in signal between the 38 pre-infiision time points 
and the 98 post-infiision time points in a fixed volume around the peak voxel with the 
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5 NAc was 1 .50%, corresponding to an effects size of = 1 .79 standard deviations. To 
achieve 90% power to detect a signal difference of this magnitude at p < 0.05 (two- 
tailed) would require N = 15 independent comparisons (Cohen, 1988). Effect sizes of 
cocaine infusion in other subcortical structures ranged from d = 0.60 to 2.14. For the 
morphine infusion study, effect sizes of morphine infusion in drug-naive subjects 
10 ranged from d = 0.71 to 1.67. We also note that our preliminary data indicates effects 
of similar magnitude can be found when comparing cocaine dependent and healthy 
control subjects with non-infusion paradigms such as the monetary reward paradigm. 
These calculations suggest that endophenotypes based on quantitative differences in 

* • * 

signal change across individuals should be distinguishable, e.g., with 15 subjects per 
15 diagnostic group. 

/ 0342] Family Based Association Studies: 

[0343] The recruitment strategy is based upon the identification of families as 

the unit of analysis. Thus we also propose family based association studies for the 
candidate gene association studies that will be performed. We have performed 
20 several family based association studies (see for example Wilk et al. 2001 ; DeStefano 
et al. 2002). 

[0344] Family based association tests (FBATs) evaluating association 

between markers and the various phenotypes will be conducted using the program 
FBAT (Laird et al. 2000). These tests are described in detail elsewhere (Rabinowitz 
25 and Laird 2000; Horvath et al. 2001). A general form of a family based association 
test statistic for family i (with nj offspring) is 

where Xy is a function of the genotype data of offspring j in family i and Ty is a 
30 function of the phenotype data of that offspring. For a bialleUc marker a score 
statistic based on Si can be defined as 
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Z = Z [S,.-E(S, )]/ Vv(S,.) 
1=7 

5 

where E(Si) and V(Si) are the mean and variance of Si under the null hypothesis of no 
linkage and N is the total number of famihes. If the coding of Xy specifies an additive 
model (i.e. Xij = the number of alleles of interest (0, 1 or 2) carried by offspring j in 
family i) and Ty is specified as 0 for unaffected and 1 for affected, then this statistic is 
10 equivalent to the TDT for genotyped parent-offspring trios (Lxmetta et al, 2000). 

[0345] When parental genotypes are available, E(Si) can be computed by 

conditioning on the observed traits and parental marker genotypes, and is based on 
Mendelian transmission probabilities (see Horvath et al. 2001 for details). This 
further justifies the collection of parental genotype information. Rabinowitz and i 

15 Laird (2000) invoke the statistical method of conditioning on sufficient statistics for 
the null hypothesis to construct a test of association when parental genotypes are not 
available. In this case the offspring genotype distribution is defined by conditioning 
on the observed traits, the partially observed parental genotypes and on the offspring 
configuration. Tables presenting the conditional probabihties when partial or no 

20 parental genotype information is available are given in the FBAT technical report 
portion of the FBAT documentation. At least two distinct offspring genotypes must 
be observed for a family to contribute to the FBAT statistic when parental genotypes 
are not available. The statistical theory of conditioning on the sufficient statistics 
results in correct p-values (type I error rate) regardless of the population admixture, 

25 patterns of missing genotypes or genetic model (Rabinowitz and Laird 2000). 

[0346] Both multiallelic and biallelic association tests can be used. For the 

biallcUc tests an additive genetic model will be assumed with Xy coded as described 
above. Coding of Xy for multiallelic tests are described elsewhere (Horvath et al. 
2001). The unknown underlying genetic model may determine which test, biallelic or 
30 multiallelic, is more powerful, hence both wiU be considered here. Two definitions of 
the trait will be employed. In the first definition Ty = the quantitative trait for 
offspring j in family i. In the second definition Ty = (quantitative trait - |a), where |li= 
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S a constant that is chosen to minimize the variance of the test statistic (Horvath et al. 
2001). For these trait definitions, a positive Z statistic indicates that the allele is 
associated with a larger value. 

[0347] Sib Pair Estimates: 

[0348] In the linkage power analysis, 900 sibling pairs are mentioned, but the 

10 analysis is for a quantitative trait, so it uses the continuous measure for fMRI as the 
trait to which we are linking. Other than the proband, sibs are not defined as 
"affected" or unaffected, merely by their fMRI measure(s), 

[0349] An exemplary 900 sib-pair is estimated based upon the following 

numbers: 



Cocaine families: 213 families, 5 members in each family: 

half (n=107) consisting of 1 parent and 4 offspring (107 probands, 321siblings), 

half (n=106) consisting of two parents 3 offspring (106 probands, 318 siblings): 

(total number of sibhngs = 639). 

Number of sibling pairs: 

107 X 6 = 642 sib pairs 

106x3 = 318 sib pairs 

960 sibling pairs total. 

Nicotine families: 152 famihes, 7 members: 

half (n=76) consisting of 1 parent and 4 children, 2 avuncular or cousin, 
half (n==76) consisting of 2 parents, 3 children, 2 avuncular or cousin* 

456 sib pairs 

228 sib pairs = 684 pairs 

For the purposes of power estimates we will add one additional sibling 
pair (e.g. one cousin pair) for each of the 152 families; 

1 52 additional pairs 
836 sibling pairs total. 

Familial Depression families: 107 families, 10 members: 

half (n=54) consisting of 1 parent and 4 children, 5 avunculsir & cousin, 

half (n=53) consisting of 2 parents, 3 children, 5 avuncular & cousin. 

456 sib pairs 

228 sib pairs = 684 

For the purposes of power estimates we will add two additional sibUng 
pairs (e.g. two cousin pairs, perhaps from two different sibships) for each 
of the 107 families: 

214 additional pairs 
898 sibling pairs total. 



15 
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5 [03501 Example (Part 3^: Detailed Description of an Exemplary Database 

[035 1 ] The text that follows summarizes a number of sahent features of the 

Brain Imaging, Genetics, and Behavioral Assessments Database (BIGBAD) including 
Database Design, Database Architecture, Data Entry Procedures, Data Transfer 
Procedures, Data Confidentiality, Accessibility and Security, Quality Control, and 
10 Database Personnel. 

[0352] (a) Database Design 

[0353] The Brain hnaging. Genetics, and Behavioral Assessments Database 

has been designed to meet the following objectives: (1) Receive and store all data 
(behavioral, MRI and molecular genetic) acquired during the project; (2) Provide an 

15 . easy to use, intuitive interface which reflects the woik- and dataflow defined by the 
project protocol; (3) Provide data entry interfaces for behavioral data entry which, to 
the extent possible, mimic the 'actual' test forms; (4) Perform immediate, automatic 
quaUty control where possible (vahdity of data entry, e.g., type, range; redundancy 
checks); (5) Provide faciUties for 'manual' quaUty control at various stages; (6) 

20 Automate data transfer (behavioral measures, MRI and molecular genetic) as much as 
possible; (7) Simplify communication between the four working cores of the 
phenotype-genotype project; (8) Serve raw and processed data to the outside world, 
under to-be-defined siccess control. 

[0354] Note that for simplicity, in this section all non-MRI and non-genetic 

25 data (i.e., clinical, neurological, cognitive, ...) is referred to as "behavioral." The 

overarching goal of data coordination is to collect all MRI, genetics, and behavioral 
data acquired for the targeted study. In the case of MRI scanning and molecular 
genetic studies, collecting data is a relatively straightforward process; in contrast, the 
behavioral data is significantly more complex, both conceptually and practically. In 
30 general, BIGBAD has been designed fi"om the assumption that, whenever possible, 
ALL behavioral data (raw data and summary scores) will be stored. 

[0355] BIGBAD has been designed in a modular fashion, making it highly 

flexible and expandable. As such, each behavioral test is implemented as a separate 
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5 database module, developed in coordination with the PI in the Clinical Phenotyping 
Working Core responsible for that instrument Many tests also required development 
of complex scoring algorithms, which were either defined or developed by the 
responsible PL The result of this common effort will be automatic real-time scoring of 
the majority of instruments at the time of data entry. Diagram 1 lists the instruments 
10 included in the behavioral test battery 

Diagram 1 : Behavioral instruments included in the Phenotype-Genotype test battery 

Commercial and/or electronic tests 

SCID -UP 

IDS 

Fagerstrom Nicotine Tolerance Questionaire 

SSACA 

SSAGA-n 

Other tests 

Full medical History, ROS and Exam 

Full neurological History, ROS, and Exam 

Handedness 

Pregnancy Test 

HIV and HepC 

LFTs, CBC, SMA-20, 

Hair toxicology 

Urine Toxicology 

WAIS-R 

Saliva Continine 

End-expiratory CO 



15 f0356J (b) Database Architecture 

[0357] The primary rationale for using an estabhshed database such as 

BIGBAD is to provide ease of communication between the working cores, ease of 
data-entry, and continuity in workflow and dataflow by closely mirroring the Project's 
logic and workflow (see the diagram for information flow in Appendix 5 that 
20 organizes the activities performed by each working core). 

[0358] The components of the system in terms of data acquisition and 

processing can include: (i) the clinical phenotyping working core and their offline 
(i.e., paper and pencil), online (i.e., computer-based), and chemistry-based meajsures; 
(ii) the MRI scaxmer used by the neiuroimaging working core and its data-analysis 
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S platforms; (iii) the quantitative anatomy working core and its data-analysis platforms; 
(iv) the neurogenetics working core; (v) the PC's or Linux-based workstations at each 
working core, that upload data to the central database; (vi) the database hardware and 
software installed and configured at the central database, allowing data entry and 
access through predefined access mechanisms; (vii) the central supercomputer and 
10 disk storage; and (viii) the data backup system (i.e., tape-farm) for the central 
database. 

[0359] The database can handle data acquisition from multiple working cores, 

provide full subject confidentiality, and manage repetitive testing/scanning of subjects 
throughout the course of the study. 

1 5 [0360] The database architecture can have a three-tier structure: 

Databzise Layer - a relational database (server side) 

Application Logic Layer -^pUcation logic controlling user access and query 

execution 

Front-end Layer -web-based graphical user interface (GUI) (for 

investigators in each working core) 

The structural three-tier organization enables applications to be distributed over many 

physical locations and computing platforms, hivestigators access the database via 

front-end interfaces (e.g., GUIs) developed to best suit their computing enviromnents. 

These interfaces can be implemented using virtually any programming language and 

20 even other databases* GUIs (for example, Microsoft Access can be used as a front-end 

to a MySQL database). At the same time, investigators can seamlessly connect to 

multiple databases using one GUI. 

[036 1 ] ( 1 ) Database software platform 

[0362] The database is developed using MySQL, an Open Source Database 

25 Management System (see http://www.mysql.com). MySQL is a database management 
system that incorporates a relational model for its databases, and supports ANSI SQL 
(standard querying language). It is very flexible and supports compatibility with other 
database management systems. MySQL also supports ODBC (Open DataBase 
Connectivity, an industry standard application programming interface (API) for 
30 transparent database access) and JDBC (a Java API for executing SQL statements). 
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5 hence making it possible to use MySQL as a back-end database to many different 
applications (e.g., MRI data processing pipeline, Microsoft Excel, Matlab, ...). 
MySQL's client/server architecture allows the development of various front-end 
interfaces with seamless connectivity to the Database servers. 

[0363] The MySQL architecture corresponds well to the requirements of the 

10 Phenotype-Genotype Study. The server (The MySQL daemon process mysqid) 

coimects investigators by creating a new server process for each investigator. 

Investigators access the MySQL database exclusively through the mysqid process. 

Thus the MySQL database server (program) focuses only on data handhng, while the 

mysqid processes take care of the investigator's connectivity and control his/her 
1 5 access privileges. 

[0364] The Graphical User Interface (GUT) to the database was developed to 

ensure data and structure flexibiUty, cross-platform independence, and transparent and 
foil Internet support. It has been implemented as a web-based GUI, written primarily 
in PHP4 (http://www.php.net). PHP is a powerful and versatile server-side scripting 
20 language, featuring an extensive programming interface to MySQL. For certain 

operations and data manipulation tasks, PHP is complemented by software developed 
in Perl, JavaScript, or Java. 

[0365] For secure and automatic data transfer from any PC/workstation in the 

working cores to the central database, a combination of Unison 
25 (http://www.cis.upenn.edu/-bcpierce/unison/) and Secure Shell (SSH, 

http://www.ssh.com) is used. Unison is a file synchronizer, which efficiently 
synchronizes the data present on the laptop with a central data repository. This 
process is run through SSH to ensure secure, encrypted data transfer. 

[0366] (2) Database Layer 

30 [0367] The core of the management system for the database is a relational 

database with thousands of fields storing neurological, psychological (behavioral), 
and medical data (including genetics data), raw and derived scores, MRI scans and 
analyzed images, and MRI header information. 
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5 [0368] The core of the database structure is the candidate profile, built around 

a study subject as a basic "data unit" A study subject is registered by an investigator 
at the chnical phenotyping working core. Some study subjects will undergo multiple 
visits to a particular working core (such as for test and retest scaiming). During such a 
visit, a distinct battery of behavioral instruments and MRI procedures will be 
10 administered, both of which are age and study objective dependent. 

[0369] (3) AppUcation Logic Layer 

[03 70] The middle tier consists of MySQL-based user management functions 

(using a special "mysql" database to manage user accounts). This enables the mysql 
daemon processes to verify user accounts at connection time and keep track of their 
15 access privileges during their work with the database. This way, the (work) load to 
verify users is removed from both Database and front-end applications. 

[0371] At the same time, PHP (server-side), Perl, and Java scripts dynamically 

develop SQL to query the Database, receive and process resulting data sets and 
present them to the front-end applications. Since the database front-end delivers 
20 completely dynamic web-content that is displayed on the investigators browsers, it's 
this appUcation layer's job to define and deliver variable and rules for displaying the 
content. 

[0372] (4) Front-end Layer 

[0373] The front-end (GUI) layer was designed to mirror the project 

25 workflow, its forms to resemble original layouts of the paper test forms, thereby 

making data entry highly intuitive. The main Menus of the GUI represent the actual 

candidate screening and data acquisition stages of the study: 1 . Candidate 

Recruitment Stage/Menu (e.g., initial recruitment of the candidate to the project); 2. 

Candidate Screening Stage/Menu (e.g., fiirther pre-visit screening of the candidate); 3. 
30 Candidate Visit Stage/Menu (e.g., candidate visit for clinical phenotyping); 

4-Approval Stage/Menu (e.g., post-visit period for evaluation of collected data and 

administered neuro-psychological instruments). 

[00374] Other Menus bring User Management, Data Management, and 
Candidate Profile Management features: 1 . Central Database Area (e.g., data 
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management features, offers real-time monitoring of the data acquisition process at all 
sites, user-defined querying and displaying of various database statistics); 2. 
Candidate Information (e.g., candidate profile management features); 3. User 
Information (e.g., user personal and contact information); 4. Administration (e.g. 
various administrative tasks, firom registering new users to changing access privileges 
for users or groups). 

[00375] The MRI and behavioral battery of instruments will be clearly 
displayed in the candidate profile menu. Data entry and evaluation status of all 
instruments will offer easy review of the status of work with each candidate enrolled 
in the project. 

/ 003 76] (c) Data Entry Protocols and Procedures 

[00377] The data entry aspect of the database has been designed along two 
principles: make data entry easy, and make it accurate. To make data entry easy, the 
online forms have been designed to resemble as much as possible the paper forms that 
researchers are used to working with. The same headers, titles and layouts as the 
paper tests will be provided online in many cases, and where they are not, clear 
instructions will be written to smooth any transitional problems. Data entry in this 
environment is extremely fast, and typically takes only a few minutes for even the 
longest measures. Many shorter measures take only moments to enter data, and 
feedback (including scores) is immediate. 

[00378] To make data entry accurate, the online forms provide several basic 
levels of quality control. They limit the entry options of nearly every field, making 
unreasonable values impossible to enter. They provide immediate feedback to the data 
entered, and allow investigators to easily check any and all of their entries. Finally, 
trained personnel will explicitly verify a randomly selected subset of the data entered 
against paper originals. 

[00379] Data entry on the commercial software integrated into the project is a 
more complex issue. Each commercial software package has its own protocols for 
data entry, but when exported information arrives at the central database, it is run 
through a standardized battery of checks. Primarily, these checks involve verification 
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of the candidate identity (does this file belong to the correct subject), and of basic 
infonnation content (does this file contain the information that it should). After these 
checks have been done, the data is subject to the same quality control as the data 
arriving via the online interface. 

[00380] Note that given Phi and Ph2 are primarily focussed on research at one 
center, the data entry system will feature double entry of data, a standard procedure 
for maintaining data quality. For Phi and Ph2, this database will perform double 
scoring of the behavioral instruments, thus providing quality control by comparing the 
summary scores only. For Ph3, this feature may be developed to the point that it can 
be utilized across multiple centers. 

[ 00381] (d) Data Transfer Protocols and Procedures 

[00382] The proposed data transfer mechanism for the phenotype-genotype 

project calls for a study workstation for each investigator in the four working cores to 
fimction as an extension of the central database, effectively constituting a data 
gateway between them. In this scenario, all acquired study data, be it 
clinicalTjehavioral, MRI, or neurogenetic measures flows firom acquisition through 
the workstation to the central database. 

[00383] Although data transfer is technically possible using currently available 
mechanisms, it should be noted that for MRI data this procedure can be cumbersome 
and require additional himian resources. Specificjally at the central database, the 
verification, QC, and format conversion of MRI data requires significant manual 
intervention. 

[00384] For data transfer purposes, this study has three primary categories of 
data: (1) clinical/behavioral data from paper-and-pencil tests and computerized tests 
(e.g., SCID-I/P); (2) structural and functional MRI data; and (3) molecular genetics 
data. 

[00385] These data types can be acquired in different ways, and may require 
slightly different treatment for storage, archiving, backup, database entry and transfer. 
In the following, the procedures around the clinicaybehavioral and MRI data are 
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described, given these are the most complex, or the largest data sets, respectively, in 
the project. 

[00386] The data transfer mechanism has all data acquired or analyzed at the 
working cores travel via a laptop/workstation (or from the MRI scanner to a 
workstation in the neuroimaging working core for offline reconstruction of images) to 
the central database. For two of the data categories, these procedures are summarized 
as follows: 

[00387] Clinical/behaviral tests: (a) one set of tests is administered to the 
subjects using standard, paper-and-pencil test forms. The data contained on these 
forms are to be entered into the Database using a data entry interface provided by 
central database, which can be accessed over the Internet. Note that data entry will not 
be limited to a single laptop/workstation; other computers at the working core, will be 
usable for data entry, (b) Another set of tests are computerized and administered using 
a laptop with a battery of computerized tests. The data generated by these instruments 
are initially stored in the internal representation of each individual software package. 
Following test administration, the test results are manually 'exported' to a format 
usable for transfer to the Database. Each laptop/workstation will be configured with 
an upload mechanism that automatically transfers the exported data to the Database. 

[00388] Structural and functional MRI data: scans are acquired at the MRI 
console, and from there pushed' to the Workstation using DICOM transfer. From the 
Workstation, they are subsequently sent to central database using a similar, encrypted 
DICOM transfer mechanism. 

[00389] (e) Data Confidentiality, Security, and Accessibility 

[00390] The database is designed with a number of features that control access 
to the database and ensure subject confidentiality. 

[00391] (f) Quality Control for Data Integrity 

[00392] Four different levels or stages of quality assurance and quality control 
have been designed into the dataflow: 
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[00393] (1) At the working core, during and after the data entry. The 

behavioral data are checked automatically for validity, type, and range as the data are 
entered in the on-screen test forms. The MRI scans are visually checked at the MR 
console. 

[00394] (2) At the working core, before the data is transferred to central 
database. For behavioral data, this includes an explicit data entry/completeness check. 
The tests are displayed in the order of administration, making it easier to monitor the 
data entry process. Once the user enters all the data for a certain instrument, he/she 
has to mark that the data entry is completed. This informs other users that the test's 
data entry is completed and disables anyone else but that user from editing the entered 
data (with the exception of the working core PI, who has the authority to access and 
modify all data of his/her working core). For MRI data, this QC stage consists of a 
qualitative evaluation of the data during pre-processing and before statistical 
evaluation, using visualization software that allows multiple simultaneous cross- 
sectional views. 

[00395] Once a test's data entry or MRI acquisition has been completed and 
checked as such, the authorized user may evaluate the instrument and mark it as 
"Completed PASS" or "Completed FAILURE". If the instrument was not 
adnainistered for some reason, it may be also checked as "Not Administered." 
Simultaneously, a record (relatjsd to the QC level) is updated, while an entry 
(comment) is inserted into the comment history table of the Database. Each time this 
is done, the QC flag table gets updated, keeping the latest entry, while the table 
comment history keeps the chronological listing of aU comments. This provides a 
complete audit trail, recording exactly what was done with the data throughout the 
course of the study. 

[00396] (3) At the central database, upon receipt of the data at it. This stage 

verifies the integrity and completeness of the received data and MRI scans, i.e. if the 
received files were correctly transmitted, whether the data is complete, and whether 
the correct acquisition parameters were used. 

[00397] (4) At the central database, following data receipt and integrity check. 
This level of Quality Control is the most comprehensive, m-depth verification of all 
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received information for a study subject. The validation at this QC level initiates the 
candidate*s "promotion" into a status of a full subject, when a study wide, unique 
Subject ID is assigned to it. For behavioral data, this involves a complete verification 
of all data against source documents (paper forms) on a random subset of candidates, 
and rapid data consistency checks of all data. For MRI data, this involves the 
qualitative and quantitative assessment of image quality 

[03 98] (g) Central Database Organizational Structure 

[0399] The central database is organized into multiple separate domains of 

activity for each of the types of data to be incorporated in it (thus approximating the 
structure for the four working cores). 
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[0400] Other embodiments of the invention are within the following claims. 
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